Probability Surveys 



Vol. 3 (2006) 413-466 
ISSN: 1549-5787 



DOI: 10.1214/154957807000000013 



On the 



constructions of the 
Brownian motion 



skew 



Antoine LejajS 

IECN, 
Campus scientifique, 
BP 239, 

54506 Vandceuvre-les- Nancy CEDEX, France 
e-mail: Ant o ine . Le j ay @ i e cn . u-nancy . frl 



Abstract: This article summarizes the various ways one may use to con- 
struct the Skew Brownian motion, and shows their connections. Recent 
applications of this process in modelling and numerical simulation moti- 
vates this survey. This article ends with a brief account of related results, 
extensions and applications of the Skew Brownian motion. 

AMS 2000 subject classifications: Primary 60J60; secondary 60H10, 
60J55. 

Keywords and phrases: skew Brownian motion, PDE with singular drift, 
PDE with transmission condition, SDE with local time, excursions of Brow- 
nian motion, scale function and speed measure, mathematical modelling, 
Monte Carlo methods. 

Received November 2005. 



1. Introduction 

The Skew Brownian motion appeared in the '70 in (44J [87] as a natural gener- 
alization of the Brownian motion: it is a process that behaves like a Brownian 
motion except that the sign of each excursion is chosen using an independent 
Bernoulli random variable of parameter p. As a sequel, if X is a Skew Brow- 
nian motion of parameter p (that is, p is the probability that the excursion is 
positive), then P[X t > 0] = p for any t > 0, while outside 0, X t behaves like a 
Brownian motion. As shown in [41] , this process is a semi- martingale which is a 
strong solution to some Stochastic Differential Equation (SDE) with local time 



where L®(X) is the local time at of A. 

At the same time, N. Portenko constructed in [BSl [70] a process whose in- 
finitesimal generator has a singular drift concentrated on some hypersurface. 
This is a way to model some permeable barrier. Indeed, he considers diffusion 
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whose infinitesimal generator is 
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+ J2<l$sN l (x) 



where N — (Ni, . . . , Nd) is the vector conormal to a smooth hyper-surface S, a, 
b and q are continuous and \q(x)\ < 1 for all x £ S. He showed that the solution 
to — Lu with u(0,x) = <p{x) is also solution to the PDE 

i (1 + q(x))N+(x) ■ Vu(t,x) = (1 - q(x))N-(x) ■ Vu(t,x) 
u(t, ■) is continuous on S, 
K u(0,x) = ip(x) 

where N + and are the inner and outer conormal to S. The condition on the 
left and right flux on the surface S is called a transmission condition or a flux 
condition. 

In dimension one with a = 1 and 6 = 0, this process is the Skew Brownian 
motion. Thus, the infinitesimal generator of this diffusion process has a rather 
natural interpretation, since it corresponds to half the Laplace operator plus a 
"generalized drift" given by a Dirac mass at with a coefficient that corresponds 
to its skewness. Besides, the effect of the drift is translated into the PDE as a 
transmission condition. It was also noted in |92j that SDEs with local time may 
be used to model diffusion processes in a media with permeable barrier. 

Since, a few works have used the Skew Brownian motion as a tool for solving 
applied problems: in astrophysics [SS], in ecology [T3], in homogenization [5DII52] . 
in geophysics [56], [57] and more recently in finance [191 [20l [2T1 122j .... 

A recent series of works [27] [56j [58j [61] shows that the properties of the Skew 
Brownian motion may be used in a systematic way to provide different schemes 
to simulate diffusion processes generated by 



with a and p have discontinuous points of first kind. All these Monte Carlo 
methods rely on a probabilistic interpretation of the transmission conditions at 
the points where the coefficients are discontinuous and appropriate changes of 
scales. 

Indeed, the diffusion process X generated by L is solution to an SDE with 
local time of type 



where v is a finite, signed measure that has a mass at the points where a or p 
are discontinuous. 




(1) 
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SDEs of type © have been studied first in [S3 ] 151 ] IS2] (see also [25] for 
some extensions) , where existence and uniqueness of a strong solution is proved 
under some rather general conditions (a shall be uniformly elliptic, bounded 
and of finite variation, b shall be bounded, and v shall be of finite mass with 
W({ ac })| < 1 for any point x £ R). 

This type of SDE generalized the usual SDEs and were studied with hypothe- 
ses on the coefficients that are as weak as possible: [HI EH 130]: ■••In addition, 
it is worth noting that the articles [531 also give an account on the "critical 
cases" between strong and weak solutions of SDEs, in the sense they deal with 
the weakest possible conditions on the coefficients to ensure the existence of 
a strong solution (see also [6] for an example of application to the theory of 
SDEs). In addition, some Dirichlet processes can be constructed as solutions of 
equations of this type [5] [35] . 

Another natural extension of the Skew Brownian motion is the Walsh Brown- 
ian motion. It is a diffusion process that moves on rays emanating from a single 
point and was introduced by J. Walsh in [87]. This diffusion process can be used 
locally as a description of a spider martingale |94[ Sect. 17, p. 103], and is a 
special case of a diffusion on a graph, since we recover a transmission condition 
at each vertex. Diffusions on graphs are particularly important in the study of 
dynamical Hamiltonian systems as shown in the pioneering work of M. Freidlin 
and M. Weber [34] . but could be useful for modelling many physical or biological 
diffusion phenomena. It is then possible to use Monte Carlo methods to solve 
that type of problems. 

One has also to note that Walsh Brownian motion was useful to provide a 
counter-example to a natural question on Brownian filtration, as B. Tsirelson 
shown it in [85] , 

The goal of this article is then to summarize the different ways to construct a 
Skew Brownian motion (using PDEs, Dirichlet forms, approximations by diffu- 
sions with smooth coefficients and by random walks, scale functions and speed 
measures, SDE, excursions theory...) and their relationships. The last section 
presents quickly some extensions of the Skew Brownian motion and their appli- 
cations to various fields. 

Notations 

Classically, we denote by R+ (resp. R_) the set of non- negative (resp. non- 
positive) real numbers, and by R* + (resp. R* ) the set R + \ {0} (resp. R_ \ {0}). 

The space of continuous functions from X to Y is denoted by C(X, Y). 

The space of square integrable function / on X is denoted by L 2 (X). 

On R d with d > 1, we denote by H 1 (R d ) the completion of the space of 
smooth functions with compact support with respect to the norm 
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A function in H 1 (M ci ) has then a square integrable generalized derivative d Xi f 
with respect to any coordinate Xi in R d . Let us recall that when d = 1, for any 
function / in H 1 (M), there exists a continuous function / such that f(x) = f(x) 
almost everywhere. 

The space H 2 (R) contains the function / in H 1 (R) such that V/ also belongs 
to H^R). 



2. Differential operators with generalized coefficients 

Let S denote the Dirac function at 0. For the first time, we are interested in 
solving the parabolic PDE for some q G [—1,1], 

^^- = iAu{t,x) + q5Vu(t,x), ^ 
it(0, x) = <p(x), 

which is equivalent to construct the semi-group associated to the operator with 
a singular first-order differential term 

L=^A + qSW. (4) 

In [551 [70], N. Portenko constructed the semi-group generated by L and showed 
it is a Feller semi-group. In fact, N. Portenko works in R m instead of R, where 
S is a Dirac mass of a surface S smooth enough. Here, we restrict ourselves to 
the simpler form of L given by ([4J, and we give a short account of the whole 
construction in Section 111.10.11 below. 



2. 1 . Perturbation of the heat kernel 



One wants to construct the semi-group (Qt)t>0i if it exists, generated by L. Let 
us denote by (Pt)t>o the semi-group generated by over R. It is well known 
that Pt has a density pit, x, y) given by the heat (or Gaussian) kernel: for any 
continuous, bounded function /, 

P t (f(x) = ( p(t, x, y)tp(y) dy with p(t, x, y) = —= exp 
Jr v 

It is natural to construct Q t as a perturbation of P t , that is 

QMx) = P t <p{x) + q f p{t-T,x,0)VQ T <p{0)dT. (5) 
Jo 

The problem here is to give a meaning to VQt<p{0), since we will see that VQt<p 
is discontinuous at 0. However, let us remark that for any continuous function 
/ and any r > 0, x h-> p(r, x, 0)/(0) is of class C 1 (R; R) and 

J-p(r,*,0)/(0)U o = 0. (6) 




A. Lejay/The skew Brownian motion 417 

So, if we inject in the right-hand side of (|5|) the value of Qt<p(x) d = u(t, x) given 
by ((5]), one gets with ([6|) that 

Q t cp(x) = u(t, x) = P t ip(x) +q f p(t-r, x, 0)VP T tp(0) dr. (7) 

Jo 

We use (O as the definition of Qt- 
2.2. Single layer potential 

Another idea to construct the semi-group is to consider that the PDE ([3]) cor- 
responds to the potential generated by a charge whose value is 

Vu(t,0+) + Vu(t,0-) 
V(t, ip) = q 

at the point 0. From standard results (See [31], [38] for example), one knows that 
the solution Q t ip(x) d = u(t,x) of such a problem is given on R^ and MI by 

Qt<p(x) = PMx) + f P (t-T,X,0)V(T,ip)dT. 

Jo 

Yet V(t, ip) also involves the value of u. Indeed, a computation similar to the 
one done previously leads also to ([7|), and thus 

V(t,ip) = VPM0). (8) 

The advantage of this formulation is that one also knows from standard results 
on potential theory (see [38] for example) that 

VQ t <p(0±)= / \/p(t,0,y)ip(y)dyTqV(r,ip) 
Js. 

which means that with {§]), 

VQM0+) = (1 - q)VPM0) and VQ^(O-) - (1 + q)VP t tp(0). (9) 

Then it is immediate that (t, x) h- > Qtf{x) is of class C 1 ' 2 on Rl x W, is 
continuous on R* x R, but x i— > VQtf{x) is discontinuous at x = 0. From ([9|), 
one gets 

aVQ t <p(0+) = (1 - a)VQ^(O-) (10) 

with 

« = — (n) 

Hence, one can give a proper meaning to ([3]), since (t, a;) i— > Q t (p(x) is a contin- 
uous solution to 

' (t,x) i-> u(t,ar) e C X ' 2 (R^ x l*;l)nC(i; x R;R), 

d"^' ^ = iAu(t,af) for any (t,x) e R* x R*, 

aVu(i,0+) = (1 -a)Vu(i,0-), ( 12 ) 
u(f,0+) = m(*,0-), 
^(Q,a:) = <p(x). 
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Proposition 1. When tp is continuous and bounded, and a G [0, 1], there exists 
a unique solution to (|12p for which the maximum principle holds. 

Proof. The existence and uniqueness of the solution of (fT2| may be proved using 
the equivalence between this problem and the same PDE written in a variational 
form: see Section l3~T1 

But we could also prove this result without using the notion of weak solutions. 
The existence of u is proved, since we have constructed such a solution with the 
help of the semi-group {Qt)t>o- 

Let us note also that the uniqueness of u follows from the maximum principle. 
If the initial condition tp is non-negative (resp. non-positive), then u is non- 
negative (resp. non-positive). Hence, if ip = 0, then u is both non-negative and 
non-positive and is then equal to 0. As (fT2")l is linear, this yields the uniqueness 
of its solution. 

If a = Q (resp. a = 1), then \12\ is equivalent to two heat equations, one 
with a Neumann boundary condition on Rl x {0} (resp. R+ x {0}), and the 
other with a lateral Dirichlet boundary on W + x {0} (resp. M*_ x {0}) which is 
specified, thanks to the continuity of u, to the value of u on Rl x {0} (resp. 
M.*^ x {0}). Hence, the maximum principle holds as it holds for both equations. 

We deal now with the case a G (0, 1). 

Assume that tp < and tp has a compact support. Indeed, if u is a solution 
of (jT2j), then u is solution to x) = §Au(t, x) both on W + x R + and R* + x 
R_, with the boundary condition ip on R + and R_, and the lateral boundary 
condition u(t, 0) on R+ x { }. The solution u for each of these equations satisfies 
u(t, x) — ► when |a;| — > +oo or t — > +oo. 

Thus, if it > inside M.* + x R (resp. R* x R) the maximum of u is attained 
at a point (t£,x£) of R^ x K + (resp. (t„ , Xq ) of R^ x R_). If x^ ^ 0, then 
it follows from Theorem 2 in [351 Sect. 1, p. 38] that u(ij,0) = u(£q ,xj). In 
other words, the maximum of u is then reached on the line W + x {0}. 

Let to be a time such that u(to,0) is maximum. As a'S7u(to,0+) = (1 — 
a)V'u(to) 0—), Vu(to, 0+) and Vu(to, 0—) have the same sign. If Vu(to> 0+) > 
(resp. Vu(£o;0— ) < 0), a Taylor development in x around shows that there 
exists some e > for which u(t ,x) > u(t o ,0) for x G (0, e) (resp. a; G (— e, 0), 
which contradicts that u(to,0) is a maximum of u. Yet Theorem 14 in [381 
Sect. 1, p. 49] asserts that Vu(to,0±) 0. Hence, u has no positive maximum 
in x 1 and thus u < 0. 

Of course, if tp > 0, then — u is solution to (|12p with — y> as an initial condition 
and u > 0. 

If </? > 0, then — < and then u > 0, since — u is solution to (fl"2")) with —tp 
as an initial condition. To conclude, let us note that if tp(x) = C for all iel, 
then u(t, x) = C. 

It remains to drop the assumption that tp has a compact support. From 0, if 
tp is continuous, \tp\ is bounded and tp{x) = on (— -R, i?), then sup^ a;) 6 R^ X R \Qtf{ 
converges to as R — > oo. By combining this result with the previous one when 
tp has a compact support, one easily gets that if tp < C (resp. tp > C), then 
the solution Q t p(x) of (jT2Jl satisfies Qttp(x) < C (resp. Qttp(x) > C) for all 
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(i,i)el + xi □ 

The next proposition follows from the previous facts. 

Proposition 2 ([69]). If \q\ < 1 (or equivalently a <E [0,1],), then (Qt)t>o is 
a Feller semi-group, where Qt is defined by (J7J) for any function if which is 
continuous and bounded on M. 



2.3. Construction of a Skew Brownian motion 

The consequences of the construction given by N. Portenko are summarized in 
the following theorem. 

Theorem 1. If \q\ < 1, then (Qt)t>o is the semi-group of a strong Markov 
process (Xt,J-t,P x ;t > 0,x G M) on a probability space (f2,.F, P), which we 
call a Skew Brownian motion of parameter a (abbreviated by SBM(a)) with 
a = (l + q)/2. This process is continuous and conservative. Besides, there exists 
a (,Ft,P)- Brownian motion B and a continuous additive functional n such that 

X t = x + B t + f}t, Vf > 0, V x -a.s.. (13) 

The additive functional r\ is of finite variation, and it variation increases only 
when the process is at 0: f Q | d?7 s | = L l{x 3 =o}l d^ s | ¥ x -a.s.. 

Of course, it remains to prove that the infinitesimal generator of X is L given 
by (@|. It is clear from ([7]) that Q t has a density q(t, x, y) given by 

q(t,x,y) =p(t,x,y) + - [ — V exp (--^ — ) dr. (14) 

t Jo V^Vt -r \ t-r t J 

In fact, a more convenient expression of q(t, x, y) will be given later in (fl"7|) . For 
that, it could be shown that for any function ip continuous and bounded, 

lim / tp(x) ( - [ q(t,x,y)(y - x) dy ] da; = #(0), (15) 



*-% V*- 

j™ / j q(t,x,y)(y - x) 2 dy j d.r = j < •(./•') d.r. (Hi) 



The equality (I16p identifies the diffusion coefficient, which is equal to 1 here, 
while (fTS"]) allows us to identify the drift term of X with q5o. This is a general- 
ization of the way of characterizing the drift term and the diffusion coefficient 
of a diffusion process, as presented for example in the book [71] . 

2. 4- An explicit construction of the semi-group 

We have constructed a semi-group (Qt)t>o such that Qtf{x) is a solution to 
the PDE (H2J). It is possible to look for a solution to (fT12|) using an explicit 
computation. 
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Let if be a continuous, bounded function. Let us remark that 
VP ^ (0) = / K J^ 6XP (^)^ )d2/ 



t->o 



= VP t £(0), 

where ^(x) = li x >o\(ip(x) — <p(— x)). Let us also remark that Ptf(x) 

when x < 0. We set ^(x) = £?(— \x\). Thus, P t 4>{x) is solution to the heat equa- 
tion on W + x W with the initial condition equal to 0. Moreover, VP t i/>(0+) = 
—WPtf(0) and VPt^(0— ) = VPt<£>(0). Thus, we are looking for a solution 
T t f{x) written under the form 

T t if(x) = P t f(x) + AP t V(x). 

where A is chosen in order that aVT t v?(0+) = (1 — a)\/T t if(Q— ). Hence, it is 
always true if 

A = 2a - 1 = q. 

This result is a direct consequence of Proposition [1] 

Proposition 3. The semi-groups (T*)t>o ond {Qt)t>a are equal. 

Hence, one gets the following expression for Q t , which is more tractable 
than (HI): 



tf(x) = / p(t,x,y)f(y)dy 



+°° 1 / 0/+H) 2 



o 

3. Using Dirichlet forms 

The idea is to get rid of the singular first-order term by transforming L into a 
symmetric operator in a proper Hilbert space. 

3.1. Weak solutions of PDE 

Here, we assume that a = (q + l)/2 ^ { 0, 1 }. The problem with the parabolic 
PDE (fT2|) is to deal with the transmission condition at given by 

aVu(t,0+) = (l-a)Vu(i,0-). (18) 
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So, an alternative way to consider (fl2j) is to look for a weak solution of the 
following PDE 



11 



'u(t,x) eC(0,r;L 2 (M))nL 2 (0,T;H 1 

^=Au(t,x), (19) 
u(0,x) = <p(x) G L 2 (M), 

where A is the divergence form operator 

A=-±-rV(a(x)V-) witho(x)= {" lfa: - + , ' . (20) 
2a(x) II — a otherwise. 

By a weak solution to (|19p . we mean a function u such that for all ip £ 
C°°([0,T];M) with ip(T,x) = 0, then, integrating formally ([19]) with respect 
to ij)(t,x)p(x) dx with p(x) = l/2a(x) and using integrations by parts, 

T / g^) 9 "^^ ^$1^ dxdt = - / <£(x)V(0,x)p(x)dx. 
o JR 02: da; 7 R 

As m is smooth on the domain where the coefficient a is smooth, it is easily 
deduced that u is smooth on (0, T] x W + and (0, T] x Ml. Let us choose ip(t, x) = 
ipi(x)ip2(t) for two functions ipi and ip2 smooth enough and such that ip2(T) — 0. 
Using an integration by parts on M+ and Ml and the freedom of choice of 
and ip2, we are led to (fl8|) . Moreover, one knows that u(t,x) is continuous on 
(0, T] x R. Thus, the weak solution of (TIT)]) is also a solution of (JT2J) and the 
converse is also true. The article 05] and the book JJSJ § III. 13, p. 224] contain 
accounts on the properties of the solution of the transmission problem. 

The domain Dom(A) of A is the set of functions / of H 1 (M) such that Af be- 
longs to L 2 (R). With this domain, A is a self-adjoint (hence closed) operator with 
respect to the scalar product (/, g)L 2 (R;a) = J R f{ x )g( x ) a ( x ) dx of L 2 (R; a). It is 
also well known that a divergence-form operator as A is the infinitesimal gener- 
ator of a Feller semi-group (P*)t>o with a density transition function (t, x, y) 1— > 
T(t,x,y) with respect to a(x)dx, that is Ptf(x) = J„T(t,x,y)f(y)a(y)dy for 
any continuous and bounded function / on M. Of course, u(t,x) — Ptip(x), and 
from the self-adjointness of A, T(t, x, y) = T(t, y, x) for all (t, x, y) G M^J_ x M 2 . 

The Proposition 0] is a direct consequence of Proposition [1] 

Proposition 4. The semi-groups {Pt)t>o and {Qt)t>o are equal. 

The next proposition follow from general results about divergence form op- 
erators: See [S3] [75] for example. 

Proposition 5. (i) As (P)t>o * s 0. Feller semi-group, it is the generator of a 
strong Markov stochastic process which is continuous and conservative (which 
is then a SBM(a) thanks to Proposition^. 
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(ii) There exists a constant 5 G (0, 1) depending only on a such that the 
function (t,x,y) <— > T(t,x,y) is S/2-Hdlder continuous with respect to t and 
5-Holder continuous with respect to (x,y). 

(iii) There exist some constants C\,C2,C^,Ca > depending only on a such 
that for all (t, x, y) G x R 2 , 

C l9 {tC 2 ,x, y) < T(t, x, y) < C 3 g(tC 4 , x, y), (21) 

where g(t, x, y) — exp(— (x — y) 2 /2t) is the heat kernel. 

The inequality (121|) is called the Aronson inequality. It allows us to deduce 
some general results about the tightness of a family of processes generated by 
divergence form operators, or a bound on the probability to leave a sphere 
centered on x prior to some time T (see for example [5U[751[75] for applications). 

Of course, the transmission condition <| 1 8|) implies that 

aV x T{t, 0-, y) = (1 - a)V B r(t, 0+, y), V(t, y) e R*+ x K. 

A note on the adjoint of A with respect to L 2 (M). When one uses 
the scalar product (/, 3)l 2 (r) = J R f( x )g( x ) dx, it is easily seen that the ad- 
joint A* of A is given by A* = -A(^-). This implies that Dom(^4*) 
{ / g L 2 (M) | f/a G Dom(A) }. Thus, any function / in the domain of A* shall 
satisfies the condition 

af(0— ) = (1 — a)/(0+) and V/ is continuous at 0. 

This also explains why the density p(t, x, y) of the SBM(a) with respect to the 
Lebesgue measure is discontinuous as a function of y, in account to the Forward 
Kolmogorov equation. Of course, this can also be deduced from Proposition[5{ii). 

3.2. Dirichlet forms 

Let (£ ,Dom(£)) be the quadratic form 

£(f,g) = ~Ja{x)Vf(x)Vg(x)dx 

for any (f,g) G Dom(£) = H 1 (R). Let (•,->l2( ,r) be the scalar product (f,g)^( a ,u) = 
Swl f( x )9( x ) a ( x ) dx. Then (£, Dom(£)) is the bilinear form associated to (A, Dom(A)) 

by 

(Af,g) LHa , m =£(f,g), for all (f,g) G Dom(A) x H 1 (M). 

It is indeed easily checked that (£, Dom(f)) is a Dirichlet form [Ml 155] , that is 
a symmetric bilinear form that is closed and contractive. 

We can now give a third construction (however with somewhat complicated 
tools that we use here in a simple case) of a SBM(a). 

Proposition 6. The Dirichlet form (£,Dom(£)) generates a continuous, Hunt 
(hence strong Markov) process, which is a SBM(a). 
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Proof. This follows from standard results in the theory of Dirichlet forms since 
(£, Dom(£)) is a regular, local Dirichlet form. Indeed, the infinitesimal generator 
of the process (£, Dom(£ )) generates is (A, Dom(A)), hence it is a SBM(a). □ 

Remark 1. In general, the process we construct using the previously invoked 
results on Dirichlet form is only defined for quasi-every starting point but not 
for any starting point. Yet in dimension one and under our assumptions on £, 
any set of zero capacity is empty. 

3.3. The Ito-Fukushima decomposition 

The theory of Dirichlet form gives us another way to prove Theorem [1] We 
denote by {X, T t , V x ; t > 0, x e R) the SBM(a) and by (6 t ) t >o its shift operator. 

A continuous additive functional (CAF) is a continuous process Y such that 
Y t+S (uj) — Y t (9 s uj) + Y s (uj). A typical CAF is given by / f(X s )ds for some 
bounded, measurable function /. 

A CAF Y is characterized by a Radon measure v, called its Revuz measure 
by the following way: If / is solution to 

(/, ¥>>Li(a,R) +£(/> <P) = [ V>(x)dv(x), V^eC^R,!), 
Jr 

then Jq + °° e~* dYt ] is a continuou^l version of /. 

Let G a be the resolvent of (£, Dom(£ )), that is the operator giving the unique 
solution u e H 1 (M) to 

a(u, v>L»(o,a) + £( u > <P) = (f; V>l 2 (q,r), V<P G H^R). 

Using the resolvent of the Dirichlet form, we see immediately that if v{ dx) = 
h(x) dx, the corresponding CAF is J Q h(X s )/a(X s ) ds. 
We say that a CAF Y is of zero quadratic variation if 

e(Y) d = sup ^ / a(x)E x [ (Y t - Y ) 2 } dx = 0. 

t^Q It J M 

Of course, a CAF is said to be locally of zero quadratic variation if there exists 
a sequence { r" }„ eN of stopping times (with respect to the minimal filtration 
generated by X satisfying the usual hypotheses) converging almost surely to 
+oo and such that N. AT n is of zero quadratic variation. 

Theorem 2 (Ito-Fukushima decomposition). Let f be a function locally in 
H 1 (R). Then there exists a local, square-integrable martingale M* with (M*) t = 
Jo l^/(-^s)| 2 ds, and a CAF locally of zero quadratic variation such that 

f(X t ) = f(X )+M( + N(. 

1 In the general case, it is only a quasi-continuous version, since a CAF may be only denned 
for any starting point except for those in a set of zero capacity. However, as noted in Remark[T] 
in dimension one and under our assumptions on E, there are no such sets. 
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Proof. This theorem follows from the results in [39]. □ 

Corollary 1. If f belongs locally to H 2 (1R) (we choose a version of f such that 
f is continuous), then 

N t = f i x s) ds + (3 fit with (if = , 



where (it)t>o is the CAF associated to the Revuz measure 6q- 
If f{x) = x, then N{ = ( ^p±i t . 

Remark 2. Of course, up to some multiplicative constant, the CAF i is the local 
time at of the process X (see Proposition VI. 45. 10 in [74]). 

Proof. Using a localization argument, we may assume that / belongs to H 2 (R), 
and we choose a version of / such that /' is continuous. By an integration by 
part, for any smooth function Lp with compact support, 

/■+oo />0 

a(x)f'(x)ip'(x)dx = - af"(x)<p(x)dx- (1 - a)f"(x)ip(x) dx 



- a/'(0+M0+) + (1 - a)/'(0-M0-). 

Hence, 

(/) ¥>)L 3 (a,R) = / (p(x)v(dx) 

JR 

with v{ dx) = a(x)(f(x) — f"(x)/2) dx — PfSo. It follows that / = G\v and that 

N(= I f(X s )ds-Y u 



JO 

where Y is a CAF associated to v. Yet the CAF associated to v is Y t = 
Io(f(X s ) - f"(X s )/2) ds - (3 fit, which yields the result. □ 

3-4- Approximations of the coefficients 

We give now two results about the approximation of the coefficients of a diver- 
gence form operator. 

Let (p, a, b) be measurable functions from K to K. We assume that there exist 
some constants A and A for which 

A < a(x) < A, A < p(x) < A and |6(a;)| < A. (22) 

Let (L,Dom(L)) be the infinitesimal generator 

Under (|22[) . (L,Dom(L)) is the infinitesimal generator of a continuous, conser- 
vative strong Markov process (See [78l [54], ... for example). 
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Proposition 7. Let (p n , a™, 6„) n6 N be a family of measurable coefficients such 
that (|22p holds uniformly and 

1 h q (R) 1 1 L 5 (R) 1 , b n L«(R) 6 

— 7 — Tj k — : — r ana 1 — (23) 

a™ n— >oo a(a;) p™ n— >oo p(x) p n a n «-+oo pa 

/or some 1 < a < oo. TTiera i/ie stochastic process X n generated by the operator 
2p™ (S ( a ™3x) "^^cEc converges in distribution to the process X generated by L. 

Proof. Under (|23|). the convergence of the resolvent and the semi-group of L n to 
that of L follows from results in [99, 100J for example. Hypothesis (|2"2"]) implies 
that one may compare the density transition functions of the processes X n and 
X uniformly with respect to the Gaussian kernel, and that these functions are 
locally Holder continuous in space and time. Thus, this is sufficient to deduce 
the weak convergence of X n to X: see (54] [76] for example. □ 

Corollary 2. Let a", p n and b n be such that (p n , a™, 6") n eN satisfies (f2"2"|) and 

(p n (x), a n (x), b n (x)) ► (p(x) , a(x) , b(x)) almost everywhere, 

n — >oo 

then the conclusion of Proposition [^] is true. 

In order to approximate the SBM by some classical SDEs, it is then possible 
to use Proposition [7| or Corollary [2] with the operator A defined in (|20| . We 
will see in the next section an application of this result to identify to different 
constructions of the SBM. We give now another application, in which the SBM 
is constructed as a renormalized limit of some SDE with drift. The following 
result has first appeared in [75] (See also Theorem [14] and Remark [15] below 
for a restatement of this result in the more general context of SDEs with local 
time). 

Proposition 8. Let X be the solution to the SDE 

X t = B t + f b(X s )ds (24) 
Jo 

where B is a Brownian motion and b G L 1 (M;R). If X" — n~ 1 X tn 2, then 
X n converges in distribution to a SBM(a) with a — e K /(l + e K ) and k — 

Proof. In dimension one, the infinitesimal generator A = 3^32 ( a g|:) ^^cE m& J 
be written 

A - — — (ae*— \ with $(x) - T 2b ^ dy 
2p dx V dxj J c a{y)p{y) 

for any £ £ R. As 6 belongs to L X (R; R), then one can set £ = —00. Assume that 
a = p = 1. The process (X, ¥ x ) generated by A is then the solution to (|24]) . The 
infinitesimal generator of X n is |A + nb(-n)4-. 
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It follows that, if $ n (x) = 2 J*^ nb(ny) dy, then $ n (x) = 2 J™^ b(y) dy and 
converges to k when a; > and to when x < 0. Thus with (j2"0|) where a(x) has 
been divided by (1 — a) and Corollary [2 it is easily checked that A™ converges 
in distribution to the Skew Brownian motion of parameter a — e K / (1 + e K ). □ 

4. Using the scale function and the speed measure 
4-1- The scale function and the speed measure 

Studying a Markov process (A, V x ) x£ r with infinitesimal generator (L, Dom(i)) 
is much simpler in dimension one, since its behaviour can generally be described 
by two unique (up to multiplicative and additive constants) strictly increasing 
functions S and V. On this description, see the books [T2l |4"4"1 174] .. .. 

The function S, defined on the state space of the process, is called the scale 
function and satisfies 

where t x = inf {t > | X t — x} is the first hitting time of the point x in the 
state space. The function S is harmonic for L, and the scale function of the 
Brownian motion is S(x) = x. A process for which S(x) = x is said to be on its 
natural scale. Of course, if S" -1 is the right-continuous inverse of S, X = S(X) 
is also a Markov process and is on its natural scale. 

The function V is a primitive of a measure m on the state space of the 
process, which is called the speed measure. Indeed, this measure characterizes 
the average exit time from any given interval for the process. If a < x < b, then 

K x[T( a .b)]= G (a , b) (x,y)m(dy) 

J a 

where Gum is the Green function 

( 2(S(x)~S(a))(S(b)-S(y)) •.f- < , l 
n n,\ - ) S(b)-S(a) 11 X ^ y> 

I S(b)-S(a) 11 X -* H 

and Tub) = m f { t > I A* ^ (a, b) }. Any Markov process on its natural scale 
may be written as B^r.), where B is a Brownian motion and T(-) is a random 
time change computed from B and m (or equivalently V). If m is absolutely 
continuous with respect to the Lebesgue measure, then the time spent by the 
process at any point has zero Lebesgue measure. 

Remark 3. The process X is the scale function S(x) = x and the integrated 
speed measure V(x) — VoS^ 1 . If the speed measure m of X has a density <p(x) 
with respect to the Lebesgue measure, then the speed measure m has a density 
^{S -1 (x)) I S' (S^ 1 (x)) with respect to the Lebesgue measure. 
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The boundary conditions are also coded by S and m, but from now, we only 
deal with processes whose state-space is R and that does not explode. 
The infinitesimal generator (L,Dom(L)) can be described by 

where 

df, > /WW 

i— (a;) = hm — ; r r ^ 

dip e-*o ip{x + e) — ip{x) 

for 95 = S or V. At a point x where ip is differentiable, j^(x) = f'{x)/(p'{x). 

The following theorem provides a way to prove the convergence of stochastic 
processes by looking at their scale functions and speed measures. It allows us 
to identify the process constructed by its speed measure and its scale function 
with a process constructed by Dirichlet forms or PDEs results. 

Theorem 3 ([36). Let (X, ¥ x ;x € R) be a conservative, continuous pro- 
cess defined by the functions (S, V), and let (X, P™;x € R)„ s n be a family 
of conservative, continuous processes defined by (S n ,V n ). If for all x G R ; 
S n (x) > S(x) and V n (x) ► V(x), then P™ converges weakly to P x in 

n — >oo n — »oo 

the space of continuous functions for any starting point x. 



4- 2. The Skew Brownian motion 

Comparing ([25]) and (l20|) suggests to construct the process (X, V x )xm with scale 
function and integrated speed measures 

„, . I oT x x if x > 0, . Tr/ . I ax if x > 0, 

S{x) = n ^-i -f Z- n and V{x) = n ^ -f Z- n (26) 
I (1 — a) x if x < 0, I (1 — a)x if x < 0. 

Remark 4. If (5 1 , V) are the scale function and the integrated speed measure 
of a diffusion process, then (kS + A, k~ 1 V + X) for k > and A G R are also 
possible scale functions and integrated speed measure. 

The question is to know if (X, P x ) x e«. is equal to the process previously 
constructed with Feller semi-groups. 

Proposition 9. The scale function S and the integrated speed measure V of 
the SBM(a) are given by (J26|). 



Proof. Let (a™, p n , &")„<=« be a family of measurable functions satisfying (f22|) . 
Let us consider 

If a™, p™ and b n are smooth enough, then it is known that A n is the infinitesimal 
generator of a continuous stochastic process (X, F™) x£ r that can be described 
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either by the functions 



5" (a;) 



o"(y) 



dyandF (z) = ^ ^ dy 



with h n (x) = 2 / s dy, 

'o P n (y)a n (y) 



or by its semi-group, which is a Feller semi-group. Hence, if almost everywhere, 
a n (x) ► a (at), p n (a;) > l/a(x) and 6™(x) > 0, then Proposition [7] 

n — >oo n — »oo n — >oo 

together with Theorem [3] prove that the process generated by A via its Feller 
semi-group and the one by its scale function and speed measure are equal in 
distribution. □ 



5. The SBM as solution of a SDE with local time 

Combining the description of the Skew Brownian motion together with the Ito- 
Tanaka formula allows us to strengthen the semi-martingale decomposition in 
Theorems [1] and [5] and to express the Skew Brownian motion as the strong 
solution of some SDE with its local time. 



5.1. The Ito-Tanaka formula 

Let / is the function from R to R which is the difference of two convex functions. 
Then#(x) = lim £ ^o !£>0 e _1 (/(x + £) - f(x)) and f'^x) = lim e _ >0 ,e<o e~ 1 (f(x + 
e) — f(x)) exists for almost every x. In addition, there exists a signed measure 
fi, called the second derivative measure, such that 

/ g(x)a( dx) = - f g'{x)f' t {x) dx. (28) 
JR Jr 

for any piecewise C 1 function with compact support on R. If / has a second 
derivative, then /" is the density of u with respect to the Lebesgue measure. 

Let X be a real-valued semi-martingale. Then there exists a process 
{L%~ (X)) t >o,x£K, called the left local time, such that for any function / as 
above, 

f(X t ) = f(x)+[ f e {X s )dX s + \ [ L*-(X)u(dx). (29) 
Jo 1 Jr 

This is the Ito-Tanaka formula. The left local time is continuous, has finite 

variation and satisfies L^(X) — J Q l{x s =x} dL x s ~ . Indeed for t > and i£l. 

the left local time L X ~(X) may be defined by 

\l>t~{X) = J o Mx s>x} dX s - (X t - x)+. (30) 

Remark 5. Here, we use a normalisation of the local time which is different from 
the one in [46] , 
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We also set for ieR, 

\L* + {X) = J !{*.<*} dX s - (X t - x)~, (31) 
which we call the right local time. For Y t — —Xt, we have that 
(Yt + x)+ = - J l {Xs <- x} dX s + \l { ~ x) -{Y) 

= -(X t -x)- = - J l {Xs<x} dX s - \L X +{X). 

It follows that L X+ (X) = -Ll~ x) ~(Y). Now, if g(x) = f(-x), then g' r (x) = 
—f'l{—x) and g' e (x) = —f r (—x). As f' r {x) differs from f' e (x) on a countable set 
[46j Problem 6.21, p. 213], the second derivative measure of g is equal to dx). 
Hence, with the Ito-Tanaka formula, 

g(Y t ) = g(-x) - J f r (-Y.) dY s - ~ J p(- dx)L x ~(Y) 
= /(*)+ f f r (X s )dX s + \ f ^{dx)L x+ {X). 

JO z JR 

As g(lt) = f(X t ), summing this expression with (f2U)) , we get 
* 1 1 '* 



f(X t ) = f(x) + J -(f' r (X s ) + ti(X s ))dX s + -J n{dx)L*{X), (32) 
where 

^ n-m+irm (33) 

is the symmetric local time. Formula (|32[) is the symmetric ltd- Tanaka formula. 
Let us note that if M is the martingale part of X, then f Q l{x 3 =x} dM s is equal 
to 0, as it follows easily from computing the expectation of the brackets of the 
integral. 

Remark 6. We have defined the left, right and symmetric local time of a dif- 
fusion process X at in view of using the Ito-Tanaka formula. An alternative 
construction is 



for the symmetric local time at 0, and 
Lt~{x) = lim- / l { x( s )e[- e ,o]} ds and L° t + (x) = lim - / 1 { x(s)e[o,e] } ds 

e— »U £ Jq e— >U £ Jq 

for the left and right local time at 0. This is a consequence of the occupation time 
formula [46, Theorem 7.1(iii)] which asserts that if / is a measurable, bounded 
function, then 

f(x)Lf(X)dx= [ f(X s )d(X) s (34) 
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for all t > 0. This formula follows easily from the Ito formula when / is smooth 
enough, and then by a density argument. 

5.2. The SDE the SBM solves 

Using the scale function and the speed measure, it is then easy to construct the 
SBM as the solution of some SDE. For this, let us consider first the SDE 

11 

Y t =y+ o{Y s ) dB s with a{y) = -l {y>Q} + l{ y <o}, (35) 

Jo a 1 — a 

where B is a one-dimensional Brownian motion and y G M. According to a result 
due to S. Nakao [63] (see also [51] [52]), this SDE has a unique strong solution. 

Thus, let B be a standard Brownian motion with its natural filtration (trans- 
formed to satisfy the standard hypotheses). The strong solution to |35|) with 
respect to B is also given by Y t = y + i?T(t) where T(t) is the random time 
change defined by 



a{B s 

The inverse r(x) of <S*(a;) is r(x) — ax if x > and (1 — a)x if x < 0. With the 
symmetric Ito-Tanaka formula (|32[) , 



r(Y t ) = r(y) + J al {Ys>0} dY s + J (1 - a)l {Yt<0} dY s + ^-l L ° t {Y), 

since J* l{y s =o} dY s = 0. Now, let us remark that 

/ al {Ys>Q }dY s + / (1 - a)l {Ys <o} dY s = / 1 {b s ^o} dS s = B t 
Jo Jo Jo 

since J Q l{ Ss=0 } d-Bs = 0. Thus, if X t = r(Y t ) = S^ 1 (B T ^ t - ) ), this proves that S 
is the scale function of X and V is the integrated speed measure of X (The func- 
tion V is identified through ([3"6"|) and Remark [3] with the help of Theorem 16.84 
in [T2] for example). 
Hence, 

X t = B t + 2 -^L° t (Y) 

for any t > 0. It remains to express L®(Y) in terms of L°(X). Again with the 
symmetric Ito-Tanaka formula, 

\X t \= [ sgn(X s )dX s +L° t (X)= [ sgn(X s )dB s + L° t (X), 
Jo Jo 

where sgn(a;) = 1 if x > 0, — 1 if x < and if x = 0. On the other hand, the 
second derivative measure of \r(x)\ is the Dirac measure Sq at 0. Hence, since 
sgnpQ) = sgn(y t ), 

\r(Y t )\ = jf* sgn(y s )r'(r s ) dF s + h° t (Y) = J* sgn(X s ) dB s + h° t (Y). 
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From the uniqueness of the decomposition of X as a semi-martingale that 
L°(X) = \L° t {Y). Thus, X is the strong solution to 

X t = x + B t + j3L° t (X) with/3 = 2a- 1, (37) 

where is the symmetric local time of X at and x = r(y). As r is one- 

to-one and since the solution to ([35]) is unique, one gets also that X is unique. 

Theorem 4 (|41j). TTie SDi? (|37p /ias a unique strong solution if and only if 
a E [0, 1]. If a strong solution exists, then it is the Skew Brownian motion of 
parameter a. 

The construction of the SDE ([3"T]) relies on the fact that S is one-to-one, S' 
is constant on R + and E* and S(R±) C R±. Then ([35]) is equivalent to ([37]) . 
We then also get that weak existence and weak uniqueness for ([3"T]) , as the SDE 
(|3"5|) has also a unique weak solution. This will be used below to deal with the 
martingale problem. 

Remark 7. In [55] , J.-F. Le Gall extended results on strong existence and unique- 
ness to SDEs of type dX t = a(X t ) dB t + J R v{ dx) dLf(X) under rather general 
conditions on the coefficient a and the measure v. In this article, the reader will 
also find another way to prove that the diffusion process generated by ([27]) also 
converges to the solution of ([37]) . 

More recently, R. Bass and Z.-Q. Chen also considered this kind of equation 
in [9] and extended some of the results of [52] . 

Section fl 1 . 81 contains a short account of this theory. 

5.3. About the left, right and symmetric local time 

As t h- ► L®(X) is a continuous additive functional that increases only when t 
belongs toZ = {s>0|A s = 0}, any continuous additive functional 1 1— > r\ t with 
the same property is almost surely equal to C n L°(X), where C n is a constant 
that depends only on ij (see Proposition VI. 45. 10 in [73]). Hence, L°(X), L°~(X) 
and L a+ (X) are all proportional. 
It follows from ^ and ^ that 

L +(X)-L° t -(X) = 2 f l {Xs=0} dX s . (38) 
Jo 

As the quadratic variation of J Q l{x s =o} dS s is equal to 0, this term is also 
equal to 0. But J Q * l {Xs=0} dL° s (X) = L° t {X). It follows that 

L° t + (X) = 2aL° t (X) and L°-(X) = 2(1 - a)L° t {X). (39) 

For alH > 0, x i— > I/^ + (X) is right-continuous with a left-limit at each point. 
Moreover, L x t ~ (X) = lim e ^o. e<o l[ x ^ + (X). Outside 0, the process X behaves 
like the Brownian motion and thus x i— > L X+ (X) is continuous on R \ { }. With 
equalities (|3"9"]) . this proves the following theorem, initially due to J. Walsh. 

Theorem 5 ([87 ). Unless for a = 1/2 (the Brownian case), the map x i— ► 
Lj + (X) is discontinuous at and continuous elsewhere. 
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6. The martingale problem 

Of course, one may wish to describe the distribution of the SBM(a) with the 
help of the martingale problem (see [46j[79] among many other books). 

Due to the presence of the local time in the SDE describing the SBM(a), this 
is not the most suitable construction. Yet we will be able to define a martingale 
problem and show that it is well posed. 

Let us denote by D(a) the set of continuous, bounded functions / on R with 
two bounded derivatives /' and /" on R* such that /"(0+) and /"(0— ) exist, 
and af'(0+) = (1 — a)f'(0— ). Hence, / is the difference of two convex functions 
and its second generalized derivative is 

1 — In 

M ( dx) = f"(x) dx + f(0-)S . 

a 

The SBM(a) is the diffusion process (X, ^F t , F x ; t > 0, x € R) which is solution 
to dX t = dB t + (3Lf(X) with j3 = 2a — 1. Using the symmetric Ito-Tanaka 
formula (El, 



./(•V, I = ./(•'■!+ / f{X s )dB s + \ ( n(dx)L*{X) 
+ §(/'(0+) + /'(0-))i?(X) 



fix) 



J o f'(X s )dB s + £f(0-)L° t (X) 



\ f f"(x)mX)dx + ^^f'(0-)L Q t (X). 



With the occupation time formula since (X) t = t, J R f"(x)L^ (X) da; — 
fo f'( x s) ds. As (3 = 2a - 1, we obtain that 

V/€D(a), f(X t ) = f(x)+ f f'(X s )dB s + \ [ f'(X s )ds 

Jo 1 Jo 

for all iel. Thus, 

M( = f(X t ) - f(x) - i J f"(X s ) ds (40) 

is a (^F^t^o-martingale under P^. 

We denote by (X t )t>o the canonical process on C(R+;R). Let us define by 
(B t )t>o the filtration B t — cr(X s ; s < t). 

Thus, let us now state the martingale problem. 

Definition 1. A probability measure Q on (C(M+; M), Bor(C(R+; R))) for which 
for all < s < t, 



/(X t )-/(X s )-i / /"(AVjdr 



'B, =0,Vf€D(a), 



is said to be a solution of the martingale problem associated to D(a). 
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Proposition 10. The distribution ¥ x of the SBM(a) is the unique solution of 
the martingale problem associated to D(a) satisfying ^ x [Xq = x] = 1 (In other 
words, the martingale problem is well posedj. 

Proof. We have already seen that the SBM(a) is solution to the martingale 
problem associated to D(a), except for the replacement of the filtration of the 
Brownian motion by (Bt)t>o- Yet, this last point may be treated as usual (see 
for example Remark 5.4.16 in [551 P- 320]). 

Let Q be a solution of the martingale problem with Q[Xo = x] = 1. 

Let us consider the scale function S of the SBM(a): S(x) — a~ l x if x > 
and S{x) = (1 - a)- x x if x < 0. Then S"{x) = for x ^ 0. The function S is 
not bounded, but with a localization argument, one gets easily that 

S(X t ) = S(x) + M t , Q-a.s., 

where M is a local martingale. Let us compute its brackets: with the Ito formula 
applied to x i— > x 2 , 

S(X t ) 2 = S{xf + 2 f S(X S ) dM s + (M) t . 
Jo 

Now, if U(x) = S 2 (x), then U also satisfies aU'(0+) = (1 - a)U'(0-) since 
U'(0+) = U'(0-) = 0. In addition, U"(x) = 2S'{x) 2 . Also with a localization 
argument, 

U(X t ) = S(x) 2 + N t + f S'(X s ) 2 ds, Q-a.s., 
Jo 

where N is a local martingale. Since 

M 2 = (S(X t ) - S(x)) 2 = 2 / S(X S ) dM s + (M) t - S{x)M t 

Jo 
t 

S'(X s ) 2 ds~S(x)N t , 

one deduces easily that (M) t = J Q S'(X S ) 2 ds. From the representation theo- 
rem (see for example Theorem 4.2 in [46j p. 170]), there exists an extension 
of the probability space (C(R+; K), Bor(C(M+; R)), Q) as well as a 
filtration (Tt)t>o and a (Q, (^ r t)t>o)-Brownian motion i? on this space such 
that M t = J Q p(s) dB s , where p is (^)t>o-adapted. Besides, Q-almost surely, 
p 2 = S'(X.) 2 . Hence, setting W t = f* sgn(p{s)S' (X s ) 2 ) dB s , one gets that, Q- 

almost surely, M t = J* S'(X S ) dW s for any t > 0. 

Since S' is constant on R + and Rl and 5(R±) C R±, one can set Y t = S(X t ) 
and then Y t = S(x) + J*S'(Y s )dW s for any t > 0, Q-almost surely. Thus, 

(Y, (!Ft)t>o, Q) is a weak solution to (|35|) . which is known to be unique. Thus 
S-^Y) = X is the SBM(a). □ 
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7. Decomposition of the excursions' measure 

Another construction is the following: Let Y be a Reflected Brownian motion, 
and let Z — { s > | Y s = } be the closure of its zeros. The Lebesgue measure 
of Z is zero, but this set cannot be ordered. However, the set R+ \ Z can 
be decomposed as a countable union U n eN Jn of intervals J n . Each interval J n 
corresponds to some excursion of Y , that is, if J„ = (£„, r„), 

Y t > for t e (£ n , r n ) and Y in = Y Tn = 0. 

Fix a € [0, 1]. At each J„, we associate a Bernoulli random variable e n which 
is independent from any other random variables (and the Reflected Brownian 
motion Y) and such that P[e„ = 1] — a and P[e n = — 1] = 1 — a. Let X be 
the process given by 

X t = e n Y t if t€J n . 

Theorem 6 (|44j). The process X is a Skew Brownian motion of parameter a. 

Proof. Obviously, the process X behaves like a Brownian motion on any time 
interval J„. By Remark 21 its scale function S and its integrated speed measure 
V are of type 

S( X ) = h" X , ifX "°' and^) = { 7 " if "-°' 
[(l- 7 )- 1 a: ifx<0 [ (1-7)2 ifa;<0 

where the constant 7 has to be specified. If T\ (resp. r_i) are the first time the 
process X reaches 1 (resp. —1), 

1 5(-l) -S(0) 

On the other hand, using the decomposition by excursions and the independence 
of the e„'s with respect to Y, 



)[n <t-i] = Y^ p o[ T i A r_i e J„; A t > on J, 
= ^ P[e n = 1 ]P [n A r_i € J n } = < 



Thus, 7 = a, and AT is the Skew Brownian motion of parameter a. □ 

Corollary 3. For all t > 0, P [A t > 0] = a. 

Let t be the right continuous inverse of the local time L®(X) of a continuous 
diffusion process X. Each jump of r corresponds to an excursion of this process, 
that is for each t such that r(i) = r(t— ), there exists an interval J„ = (^ n ,r n ) 
such that = T (t—) an d r„ = r(t). Let W be the set of excursions, that is the 
set of continuous functions / from R + to R such that /(0) = 0, f(t) = for all 
t > C for some < > and either /(*) > for t e (0, C) or f(t) < for t G (0, C). 
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For any n G N, f(t) = X((t — r n ) A (£ n — r n )) is some element of U, with the life- 
time C = V-n — r n . Accordingly, for each t € J where J = { t > ] r(t) 7^ r(t— ) }, 
one may associate a point /(t) in £/, corresponding to some excursion of X. A 
striking result due to K. Ito is that the process (t, f(t))t>o is a homogeneous 
Poisson point process with values on R + x U and intensity measure dt x P. The 
measure P is a cr-finite measure, but which is infinite. This measure P, called 
the excursions' measures, fully characterizes the diffusion process. 

For a measurable subset T of U such that P[T] is finite, F[T] denotes the 
average number of excursions in T per unit of local time. 

We denote by U + (resp. U~) the set of positive (resp. negative) excursions, 
i.e. the excursions / for which f(t) > (resp. f(t) < 0) on (0, C)- 

On the excursions theory, see pT| [74]... 

The following proposition is a direct result from the construction of the Skew 
Brownian motion given by the theory of excursions (see also 

Proposition 11. For the Skew Brownian motion of parameter a, the measure P 
may be written 

P = aP+ + (1 - a)F~, (41) 

where P + (resp. P~ ) is the excursions' measure of the Reflected Brownian mo- 
tion on M. + (resp. on that is ( | Wt | )t>o ^ (resp. (—\Wt\)t>o), where W is a 
Brownian motion. 

Remark 8. In [89] (see also [43]), S. Watanabe shows how to construct a diffusion 
given its excursions' measures, which provides another way to construct the SBM 
from Proposition mi 

Remark 9. This proposition has to be connected with the relations between the 
left, right and symmetric local times in ([39]) . 

Remark 10. For a diffusion Y with a scale function S Y , The function l/S Y (x) 
is also the average number of excursions reaching the point x per unit of local 
time. Thus, this proposition is also coherent with the fact that if S B (x) — x is 
the scale function of the Brownian motion, the average number of excursions 
of the Skew Brownian motion reaching x > (resp. x < 0) per unit of local 
time is equal to a/x (resp. (1 — a)/x), and thus the scale function of the Skew 
Brownian motion is S(x) given by ([26]) . 

Remark 11. The article [55] studies the link between the coefficients in a de- 
composition of type ([4T]) for a general diffusion with coefficients discontinuous 
at 0, that is related to the derivatives of some Green functions at 0. We recover 
the decomposition (|4"T]) with L = 2a ( a ) ^ { a ( x )'S^) with a(x) = a if x > and 
a(x) — 1 — a if x < 0. 

Another construction of the density transition function. The reflection 
principle for the Brownian motion together with the decomposition of excursions 
allows us to give another construction of density transition function of the Skew 
Brownian motion, following a result of J. Walsh [57]. Let r be the first time the 
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Skew Brownian motion hit 0, and B be a Brownian motion. Then 

F x [X t £ dy]=F x [X t £ dy;T<t}+¥ x [X t e dy;r>t]. 

If (x,y) belongs to R* + x R_ or to K* x R + , by the continuity of the path, 
P x [X t e dy;r > t] = 0. If (x, y) belongs to R* + x R + or to R*_ xL, then 

P, [ X t £ dy; r > t } = F x [ B t € dy; r >t } = (p(t, x,y) - p(t, x,-y)) dy, 

where p(t,x,y) is the transition density function of the Brownian motion. On 
the other hand, by the reflection principle and the construction of the Skew 
Brownian motion, 

P x [B t e-dy] ifx>0, y>0, 

d-aWAB^dy] ifx>0, y <0, 

^aF x [B t e dy] if x < 0, y > 0, 

Xl-a)F x [B t G -dy] ifx<0, y<0. 

Hence, after a short computation, this gives us Formula (|17j) . 
8. Approximation by random walks (I) 

As one can expect, a Skew Brownian motion may be approximated by the 
following random walk: Let (Sk)k>o be the random walk starting from with 
probability transition 

[Sfc+i =i+l\S k = i] = F[S k+1 =i-l\S k = i] = iifi^0, 
[Sk+i = 1 I Sk = 0] = a, 
[Sfc_i=-l|5 fc = 0]=l-a. 

We set for alH > and any integer n, 

1 n 2 t- Ln 2 tj 

-^■i — —b[n 2 t] H (.Jl+Ln^tJ -JL" 2 <J/- 

The following theorem may be found in [41] , and then in |52j and in [17] in a 
more general setting. P. Etore used this result in |27j to construct a numerical 
scheme and compute its speed of convergence (see Section lll.9.3[) . We give 
another proof of this theorem in Section [9] by constructing the random walk 
(S^)keN from the trajectories of the SBM(a). 

Theorem 7 f[41j). The sequence (X n ) ne jq converges in distribution in the space 
of continuous functions to the SBM(a) X. 

Proof. Let us prove first the convergence of the marginals of X n . For this, we 
remark that 

ns k =t] -- 
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Hence, 

E[exp(iAX { n )] = cm// 1 (A) + (1 - a)ijS\-\) with V(A) = E[exp(iA|X™|) ]. 

Using the results of Section [7J one has also that 

E[exp(iAX t )] = cr^(A) + (1 - a)<ip(-X) with V(A) = E[exp(iA|X t |) ]. 

The results follows from the convergence of the normalized reflected random 
walk | AT™ | to the reflected Brownian motion \X\ given by the Donsker theorem. 

The proof of the convergence finite-dimensional distributions of the X n, s is 
similar and uses the same kind of computations as in (|42|) . Yet the computations 
become heavy so that we skip it: in the next Section, we will see how to construct 
our random walk (S^)ken from a trajectory of X, and then how to prove the 
convergence of AT™ (which is then defined on the probability space of A") in 
probability to X t . The convergence of the finite-dimensional distributions is 
then immediate. 

The tightness of (Af") n6 jq follows from the Kolmogorov criteria [46j Theo- 
rem 2.8, p. 53]. For this, let us note that 

|A7-A7| 2 < A(n 2 i- \n 2 t\) 2 {S 1+ ^ t \ - S [nH] ) 2 

Let us set = Sk+i — Then all the are independent with variance 1, 
and E[£fc ] is equal to 2a — 1 if Sk = and to if Sk 0. Thus for any integer 
P > 1 > o, 

p-i /p-i \ 2 

E[ (S p - S q f } = Var(6) ^ E - 

k=q \k=g J 

< V - q + (2a - l)(p - qf < 2a(p - q) 2 . 
Hence, it follows from that for any s,t > 0, 

E[ | AT™ - Xl l \ 2 } < 3max{2a, 1}(* - s) 2 . 
One deduces that for any 7 < 1/2, there exists a random variable K™ such that 
sup \X? - X 7 s l \ < K2(t - sy and su P E[lC ] < +00. 

0<s<t<T n£N 

This implies the tightness of (AT™) n£ N. □ 

Remark 12. It follows from this proof that the trajectories of the SBM(a) are 
7-Holder continuous for any 7 < 1/2. Proposition 2.1 in [14] proves the following 
fact: almost surely, the modulus of continuity 

w(X,6)= sup \X t -X s \ 

s,te[0,T], \t-s\<5 

of the strong solution X on [0, T] to (|3T[) is smaller than 2w(B, S) for any 5 > 0. 
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In [17], it is proved that (X™, Y")t>o converges in distribution to {Bf, 
Jo f{Bs ) dB" )t>o, where (Jf") t >o and (Y t ") t > are the linear interpolations of 
= iV^r'Sk/n and = E- =1 f{Xl_ x)/n ){X- /n -X^_ l)jn ), and / is a 
measurable, bounded function. This allows us to prove some results related to 
the "horizontal-vertical" random walk. 

Remark 13. In 41j, it is also noted that if given Si — 0, the distribution of Si+i 
has the distribution of an integrable random variable Z with values in Z, then 
rr l Sy nH!i converges to a SBM(a) with a = E[ Z+ ]/E[ \Z\ }. 



9. Approximation by random walks (II) 

We give another way to prove Theorem [7] that starts from a trajectory X(ui) of 
the SBM(a). 

For the sake of simplicity, let us assume that X = 0. Fix some integer n, and 
let us construct recursively from X the following sequence of stopping times: 
tR = and 



'fe+ 



t = inf < t> T t 



X t = — for some integer £ ^ £' will X r n = — 

n k n 



Of course, as the trajectories of the SBM(a) are continuous, \X T n j — X T n\ — 

n~ x . In other words, the sequence {r]})k=o,i,... records the successive passage 
times of X(lo) on the grid { k/n \ k e N }. Finally, set 5^ = X T ^ and define X n 

by 

1 f — kSt 

X? = -SI + ——(S% +1 - S£) when t € [kSt, (k + l)6t] 
where St — 1/n 2 . 

With the scale function given in ([26]) . P[X T ™ i = (£ + l)/n \ X T n ] is equal 
to a if £ = and to 1/2 otherwise. With the strong Markov property of X, this 
means that the random walk (5^)fcgN is equal in distribution to the random 
walk previously constructed in Section [5] 

This theorem is a "specialization" to the SBM(a) of Theorem 4.1 in [52"] . 

Theorem 8. The process X n constructed from X as above converges uniformly 
on [0, T] in probability to X with respect to P. 

Proof. We prove first the convergence in probability of X" to X t for any t S 
[0, T]. Let us note first that since the reflected excursion of the SBM(a) are the 
same as the one of the Brownian motion, then for k = 0, 1, . . . , the increments 
tJ} +1 — tJ} are independent and have the same distribution (whatever the value 
of a). In particular, E[r^ +1 — ] = 1/n 2 = St (this explain our choice of St). 
The key relation is 

yn 1 on y 

A L« 2 *J/« 2 _ n l nHi ~ "" 2 *J' 

Hence, 

P[|X t n -Xf| > C] <F[\X?-X?n | > C/2] +F[\X t -X T n \>C/2\. 

[n 2 tj Y n *J 
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In the proof of Theorem [Jj the tightness of (X n )„ S N has been established by 
showing it satisfies the Kolmogorov criteria. Hence, if osc(X n : 5) is the modulus 
of continuity of X n , then for any 5 > and any 7 < 1/2, 

P[|X t " - 2 I > C] < F[osc(X", J) > C; |Tfr tJ - *l < <*] 
+ P[osc(X™,<5) > C;|if n2tJ - 1| ><J] 

<^E[^]+P[|r5, atJ -t|>H (44) 

where K" is the random 7-Holder constant of X n which is known to satisfy 
sup ne ^E[K™ ] < +00. One may replace X n by X in (|44"|) . and the result follows 
from the convergence in probability of r" n2t j to t. By a scaling argument, (r^ +1 — 

Tfc )fc=o,...,|ra 2 t| is equal in distribution to (n er*; )&=(),... i n 2tj , where the cr^'s are 
independent, identically distributed with the distribution of the first exit time 
from [—1,1] of the Brownian motion starting from 0. The Ofc have a mean equal 
to 1 and a finite variance, so that the desired result is a consequence of the weak 
law of large numbers. 

Now, to prove the uniform convergence in probability of X n to X, it is suf- 
ficient to note that for = ti < <2 < • • ■ < ti — T with t. L+ i — t.- L < S for 
i = l,...,t-l, 

sup \X?-X t \ = sup sup \X»-X t \ 
te[o,r] i=i,...,£-ite[ti,ti+i] 

< sup ( sup \X?-X£\ + \X? i -X u \ + sup \X tt -X?\) 
i=i,...,i-i \te[t it t i+ i] te[u,U+i] J 

<osc(X n ,6) + osc(X,5)+ sup \X£-X U \. 

i=l,...,£-l 

Hence, one deduces that linin^oo P[sup tg r T i |X t n — X t \ > C] = from the 
stochastic equicontinuity of (JT n )„ S N and of X by choosing 5 (and then the 
points {ti}f =0 ) small enough, and afterwards using the convergence in probabil- 
ity of the X t ™ to X ti for i = 1, . . . ,£. □ 

In order to set up a simulation scheme, the rate of convergence is computed 
in [2JJ HE] ! in the more general context of diffusion processes with discontinuous 
coefficients (besides, using a proper time increment, we are not bound to use a 
random walk on a regular grid as shown in [28] . See also points D. and E. in 
Section HH3}. 

Proposition 12 (|271 [28] ). For any 7 < 1/2 and any T > 0, there exists a 
constant C~ i ^t such that sup tg [ r] ^[ \%t — \ ] — C^^n 1 . 

10. A "follow the leader" construction 

We present here another construction based on techniques coming from the 
theory of continuous multi-armed bandits (See [60] for example). Indeed, we 
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present here in a simpler case the proof given in [2] that deals to a variably 
skewed Brownian motion (See Section 111. 6[) . Some of the results given in [T7] 
are also related to this construction. 

Let us set a € (0, 1), /3 = 2a - 1 € (-1, 1) and 7 = (1 + /3)/(l - (3) > 0. 

The idea is to construct a process (Z^ , Zf,U^)t>o with values in (K+) 2 x 
(M + ) 2 x K + from which one can get the positive excursions of the SBM(a) from 
Z 1 , its negative excursions from Z 2 and its local time from U. 

Indeed, Z = (Z l , Z 2 ) moves either horizontally or vertically when away from 
the line T = { (x, y) \ y — -fx }. The horizontal (resp. vertical) displacements of 
Z away from (U^jUl) correspond to a positive (resp. negative) excursions of 
the SBM(a) arising when its local times takes the value L t — (1 + "f)Ul . So, 
we can read "graphically" the signs and the heights of the excursions on the 
local time scale: See Figure [TJ It is then possible to reconstruct a SBM(a) from 
(Z\Z 2 ). 

{uliul) 





Fig 1. A representation of the excursions of the SBM(a) 



The heuristic idea is to construct Z 1 and Z 2 as two time-changed indepen- 
dent Brownian motion (-B-r 1 (t))*>o an( ^ (^^(t))*^ wnere either T x (i) or T 2 (t) 
increases at rate 1 while the other remains constant. The times at which the 
switches between T 1 and T 2 occur are obtained by comparing the supremum, 
as a function of time, of 7-B 1 with the ones of B 2 

We now turn to the rigorous construction. 

Let us define a subset D of (R+) 2 by 



D={ (si,s 2 )G 



7 sup B l {ui) > sup B 2 {u2) 

«l<si v. 2 <s 2 

whose closure D has the following properties: (a) R + x { } C D; (b) {s £ 
D} e T s , where (^^^(s^em 2 is the multi-parameter filtration defined by 
F{.,t) = H V T 2 - (c) If ( Sl , 8 3 ) G D, then 

{ {ui,u 2 ) S R 2 + 1 111 > si, < u 2 < s 2 } C D. 

It follows from these properties that D is described by the part of its boundary 
that intersects (M^_) 2 : See Figure [2j In addition, there exists a unique process 
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T(t) = (T x (i),T 2 (i)) — called a strategy — such that [88] (a) T(t) = (0,0) and 
T 1 , T 2 are non-decreasing; (b) T 1 ^) + T 2 (t) = t for t > 0; (c) for all s'i, s 2 > 0, 
the set { (si, S2) | < si, T 2 (i) < s 2 } is ^"( SliS2 ) -measurable; (d) the graph 

of (T(t)) t >o parametrizes the boundary of D. 




Fig 2 . Construction of the time changes 



Indeed, T x (i) increases at rate 1 when B^ 2 ^ > 7-B^w t ) and T 2 (i) increases 
at rate 1 when B%, 2 m < jB^uy m °ther words, the strategy "follows the 
leader" between jB 1 and B 2 . 

Let us now consider a maximal time interval (£,£') with t' > t such that T 1 
is strictly increasing on it, which means that T 2 is constant on (t,t'). We first 
note that as (T 1 (t) , T 2 (t)) belongs to the boundary of D, 

for any u G [t,t'], 7 sup B\ = sup B 2 2 . 

Sl e[o,ri(«)] S2 e[o,r=(t)] 

As (i,t') is chosen to be maximal, for all u G [i, i'], su P s1 g[o,t 1 (u)] -^li — ^t 1 ^') = 
Bpjjj, which means that (-Bjn(«))«G[t,t'] performs an excursion below the level 
B^ 1(t) on the [t,t'}. As (T 1 ^), T 2 (t)) belongs to the boundary of D, there 
exists no time t < T 2 (t) such that sup sg [Qt\B 2 = B^ ^ . This means that 
£?T 2 (t) = sup t6 [ 0T 2( t )] B 2 . Hence, at time time t, B^,^ has a local maximum 
and B^, 2 ^ reaches its maximum. 

We now set for r > 0, Z\. — B l Ti ^ for i — 1,2 and U} = sup sg j s j i? x , 
{J 2 = 7^. From the previous considerations, 

7 C/ 2 = 7 E# = 1 Z] = 1 Z], - Z 2 = Z 2 = f/ 2 = c/ 2 . 

Thus, on the time interval [t, t 1 ] above, Z 1 performs an excursion below the level 
XJ\ and then Z moves horizontally away from T, with Z t = Zy = Ut G T. 
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A similar study could be done for T 2 , and we see now how to define our 
candidate to be a SBM: 

X t = {ill} - Z 2 ) - (Ul - Z\) and L t = (1 + 7 )[/i. (45) 

Proposition 13. The process X is a SBM(a) and L is its local time. 

Proof. Let (Gt)t>o be the filtration defined by Q t — J~T(t)- Let us remark first 
that B t = -Zf + Z\ has for quadratic variation (B) t = T x (i) + T 2 {t) = t from 
the very definition of the strategy and is (£ t ) t >o-adapted. By the Paul Levy's 
theorem, B is a Brownian motion and then X t = B t + f3L t with (3 = 2a — 1. 
The process |X| is equal to the distance of Z from T and then, by construction, 

\X t \=B t + L t 

with B f = — B^x^ — B^, 2 ^ t y Again with computing the quadratic variation of 
this latter process, B is a (C/ t ) t >o-Brownian motion. Besides, 

-B t = Z}+Z 2 = [ sgn(X s )dB s = f sgn(X s )dX s 
Jo Jo 

where sgn(a;) = 1 if x > 0, sgn(ir) = — 1 if x < and sgn(x) = if x = 0, 
as J* l{z,^Y}d£ s = if Xt 7^ 0. Thus, with the Ito-Tanaka formula, L t is the 
symmetric local time of X. Finally, from (|45|) . X t = B t + (7 — l)/(7 + l)Lt, and 
the skewness parameter a is identified from value of 7. □ 

11. Applications and extensions 

11.1. General constructions of "Skew" diffusion processes 

The general theory of one-dimensional diffusion processes gives all the tools to 
construct rather easily a "skew" diffusion process. There are basically two ways 
to do so. 

First, consider diffusion processes living on [0, +00), or symmetric processes, 
and the sign of each excursion is changed with independent Bernoulli random 
variables. For this, one may use the description of a strong Markov process in 
term of a minimal process {i.e., a process that is killed when it reaches 0) and 
an entrance law (See [11] for a whole account on this theory) . See 86J for related 
results and constructions. 

In [H [5] a Skew Bessel process is constructed this way. 

Indeed, the symmetry of the process plays an important role in the previous 
construction, so here is another one. Let X be a process with a scale function 
S and a speed measure m, xq in the state space of X and 7 > let us define a 
new process A 7 by the scale measure S 1 and speed measure m 7 defined by 

S^(x) = h S( - x) iiX - X0 ^ndm^dx) = h~ lm{dx) iiX ~ X0 > 
\S(x) if x < xq, \m(dx) if x < xq 
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if S is such that S(xq) — 0. As the behavior of a process around a point x 
is specified locally by its scale function and speed measure and the process 
generated by (S,m) is the same as the process generated by (XS, A _1 m) for all 
A > 0, the process X 1 behaves like X except at Xq. The "skewness" of this 
process may be appreciated by the result stated in Remark [TOl 

The article |21] shows how to construct an asymmetric Bessel process by this 
method and provides applications to the computation of Asian options with 
discontinuous volatility. In [55], various processes are constructed by this way, 
in order to model some prices whose volatility suddenly increases or decreases 
when the price goes above or below a threshold. 

11.2. Walsh's Brownian motion 

Initially introduced by J. Walsh in [87], the Walsh's Brownian motion is a 
diffusion process on a set of n rays in M 2 emanating from 0. To each ray li 
is associated a weight oti corresponding heuristically to the probability for the 
process to go in this ray. On each ray, the process behaves like a Brownian 
motion. 

Of course, due to the irregularities of the trajectories of the Brownian mo- 
tion, this description is a non-sense, but this process may be described by its 
excursions' measure 

n n 

p = with X! = 1 and a * e (°> 

i=i i=i 

where P l is the excursion measure of the Reflected Brownian motion on the ray 
Ii. On this process, see also 0]. 

The Skew Brownian motion corresponds to the case n = 2. The spider mar- 
tingale generalizes the notion of Walsh's Brownian motion. This object has 
given rise to an abundant literature on Brownian nitrations, especially by giv- 
ing a negative answer to the question "if a Brownian motion is adapted to some 
filtration, is this filtration generated by a Brownian motion?" (See [HIES], [SH 
Sect. 17, p. 103] and the cited works within), but this goes far beyond the scope 
of this article. 

1 1.3. A generalized arcsine law 

For the Walsh's Brownian motion of Section fll.2[ let A l (t) be the occupation 
time up to time t of the ray 7 i; that is A l (t) = J Q 1{ x(s)eii } ds- As for the 

Brownian motion, A l (t) d =' tA l (l) for any t > 0. The arcsine law for the 
Brownian motion may be extended to this case. 
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Theorem 9 (|5J). Let T\, . . . ,T n be independent, stable random variables of 
index 1/2. Then 



E,= l a 3 T 3 



i— l,...,n 



More is proved in [5], and Skew Bessel processes are also considered. 

Assume that n = 2 and that T±, T 2 are two independent copies of a stable 
random variable with exponent 7 € (0,1) (i.e., E[exp(— (Tj)] = exp(— £ 7 ) for 
any £ > 0). Define by £ 7iQ , the random variables 



a 1 hT 1 + (l-a) 1 hT 2 ' 
The distribution function of £1/2, a is 



P[CiA« <x} = ~ arcsin y ^ + ^ _ a))a(i _ ^ . 

For a = 1/2, this is the arcsine distribution. 

In Theorem [9l the occupation time of M. + or M_ at time 1 of the Skew 
Brownian motion is a random variable of type £1/2, «■ 

If the Skew Brownian motion is replaced by a Skew Bessel process of dimen- 
sion 2 — 2v, then the occupation time at time 1 is distributed as one of the £„ ;Q 
[H [SD]. In [5U], it is also proved that for a general one-dimensional diffusion, 
the possible limit of A^+(t)/t, where A^+(t) is the occupation time of R + , is a 
random variable of type ~ . 

In the case of a Brownian motion on a graph (See below Section ril.9.2|) . 
analytical considerations allows us also to compute the Laplace transforms of 
the occupation time on each of the edges of the graph [23 . 



11.4- The Skew Brownian motion as a flow 

As the SBM is a strong solution to the SDE, 

XI = z + B t +(3L° t (X z ) (46) 

it is natural to study the dependence of the SBM(a) with respect to its starting 
point. 

Proposition 14 ([3]). For x, y £ M, let X x and X y be the two strong solutions 
of (|46p with the same parameter (3 G [—1, 1]\{0} and the same Brownian motion 
B. Then almost surely, there exists a finite time t such that X£ = X\. 

From this coupling property, we deduce the following: for /? € (— 1, 1) \ {0} 
and 1 / t/ £ 1, there exists almost surely some r > such that x + (3L T (X X ) = 
y + [3L T (Xy). 
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If Kf (X x ) = x + f3L t (X x ), the uniqueness of the strong solution to (|47| 
implies that K X S (X X ) = Ry{X x ) for s>t. 

From this property, K. Burdzy and Z.-Q. Chen give in [13j several results 
related to the distributions and behavior of the processes x G K i— ► L®(X X ) 
with t = inf { t > | L t °(X°) = 1 } and (3 € Q n (-1, 1) ^ (3L° t (X°). 

It is also possible to extend (|46|) to i € R by defining the Brownian motion 
(-B t ) te N as £? t = lR + (t)W t + 1r* (t)W^_ t for two independent Brownian motions 
W and W'. Let us note that t i— > L^A 21 ) is a non-decreasing function on M. 
If X 31 is the solution to (|4^|) . then = X x t is solution to ([4^)1 with £? t replaced 
by B^ t and /3 by — f3. Then iff (F 1 ) also satisfies the coupling property on R + . 
This implies that almost surely, there exists some time r < such that for t < r, 
Kf{X x ) = if^A^) and iff (I s ) ^ A7(A^) for t > r. 

This phenomena is studied in [131 1 14] , where the authors consider the follow- 
ing generalization of (14^1) , which is for s,x S K the solution [X s ' x , £ s ' x ) to 



Among many results, it is shown in [2] that (|47[) may have for some (s, x) € K 2 
three distinct strong solutions. The notion of "lenses" allows us to characterize 
the points (s,x) for which the decoupling/coupling of K X (X X ) takes place. 

It is then conjectured that, although there exists almost surely a strong solu- 
tion to Af = x + Bt + (3L®(X X ) simultaneously for all x £ Q, there is no strong 
uniqueness for a solution defined simultaneously for all x € R. We refer the read- 
ers to the articles cited above for a complete account on these developments, 
results and proofs. 

11.5. Comparison principle 

A natural question is to know whether or not a comparison principle holds for 
the SBM. We give first two propositions. The first one concerns the case where 
the skewness parameters (1 + /3)/2 are the same, but gives a strong statement 
on the solutions to ([47f . 

Proposition 15 ( |13| Proposition 1.8]). For all (s, s' , x, x') <G Q 4 simultaneously 
with a fixed [3 e [-1,1] \ {0}, either X*' x < xf' x ', felor X^ x > xf' x ', 
i€R almost surely, where X s ' x is solution to (j4T[> . 

In particular, if x < y, then Xf < Xf t > 0, where X x is the strong solution 
to (|4"c?|) . The second proposition concerns the case of different parameters. 

Proposition 16 ([551 Theorem 3. 4], [96 ). Let X 1 and X 2 be two strong solu- 
tions to (|46p respectively with starting x 1 and x 2 and parameters Pi and /?2- If 
x 1 < x 1 and /3i < (3 2 , then Af < X? for allt>0. 




(47) 



The distance between two solutions may then be computed. 
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Theorem 10 ( 96J). Let X 1 and X 2 be two strong solutions to (|46|) with either: 

(a) The same starting point x and two different parameters ft\ and P2 in (—1,1). 

(b) Different starting points x 1 and x 2 and the same parameter (3 G (—1, 1)\{0}. 
Then, the L 1 distance between X 1 and X 2 is given by either 

(a) E[ \X] - X 2 \ } < - j3 2 \ [ -^=I(t, x) with I(t, x) = exp ( ds, 

Jo V2tts V 2s/ 

(E[ \X} - X 2 \ }<\ x i-x 2 \ + |/3| • \I(t, x') - I(t, x 2 % 

l E [ \ L °t( Xl ) - L t( X ")\ ]<wk 1 -^l + x 1 ) - I(t, x 2 )\. 



11.6. Generalization of the Skew Brownian motion (I) 

A possible generalization of the Skew Brownian motion is to use a variable 
"skewness coefficients" , which leads to an SDEs with local time of type 

X t = x + B t + f3{L° t (X)), (48) 

where /3 is a differentiable function on R + with f3'(x) € (—1, 1) for all x € R+. 
This SDE has been studied in [2]. 

Theorem 11 ([!])• There exists a unique strong solution to (|48|) . 

The proof of the theorem relies on an extension of the construction given in 
Section 1101 The variably skew Brownian motion has been used recently for the 
Skorokhod embedding problem [18]. 

11.7. Generalization of the Skew Brownian motion (II) 

Another natural generalization is the solution of the SDE 

X t = x+ f a(X s )dB s + f b(X s )ds+ f a(s)dL° s (X), (49) 
Jo Jo Jo 

which is motivated by some homogcnization — i.e., change of scale — re- 
sults [92] . 

Theorem 12 ([91]). If a is bounded below by a positive constant, bounded and 
Lipschitz continuous, b is bounded and a is continuous with \a\ < 1, then the 
solution of (|49[) is pathwise unique. 

Remark 14. In [51] , S. Weinryb uses a normalization of the local time different 
from ours. 

11.8. Generalization of the Skew Brownian motion (III) 

The SDE solved by the Skew Brownian motion may be generalized as a SDE of 
type 

X t = X + [ a(X s )dB s + [ v(dx)Lf(X), (50) 
Jo Jr 
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where v is a signed measure. Of course, the case a — 1 and v — (2a — l)8o 
corresponds to the SBM(a). 

Let BV be the space of functions / : K — ► M of bounded variation such that 
i) / is right continuous, ii) There exists e > such that f(x) > e. 

Let M be the space of all signed measures vonl such that < 1 for 

all i£l 

Theorem 13 ( [21 [2H [STJ H2] ) • There exists a unique strong solution to (|50[) 
when a € BV and v £ M. 

Of course, the proof relies on an argument similar to the one given in Sec- 
tion 15.21 on the SDE that the SBM(a) solves. Here, the idea is to construct a 
harmonic function by the following way [52J, Lemma 2.1]: for v G M, there exists 
a unique (up to a multiplicative constant) such that 

f(dx) + (f(x)+f(x-))v(dx)=0 

and, with the additional assumption that f(x) — > 1 when x — > — oo, 

/„(*) d ^ /(*) = exp(-2^((-oo,x])) [] YTWTY (51) 



As / is in BV, f'(dx) is the bounded measure associated to / as in 
Moreover, v c is the continuous part of v. In some sense, the measure v corre- 
sponds to the logarithmic derivative of the function f(x). Then, we set F v (x) = 
Jo f»( x ) ^ follows that X is solution to ([50]) if and only if Y = F V (X) is 
solution to 

Y t = F V {X Q ) + f (a/,)(F~ 1 (Y;))dB s . 



The existence of a weak solutions follows from the martingale problem [791 146] 
while the pathwise uniqueness follows from a result due to S. Nakao [S3]. The 
Yamada-Watanabe theorem allows us to conclude (See [IB] for example). 

For the SBM(a), F v (x) corresponds, up to a multiplicative function, to the 
scale function S(x) given in ([2"o]) . 

In [52] , J.-F. Le Gall also provides the following convergence theorem. 

Theorem 14 ( 52, Theorem 3.1]). Let (^ n )neN be a sequence of measures in 
M and (<7 ra )neN be a sequence of functions in BV such that for some constants 
M and e, i) sup„ eN |^„(K)| < M, ii) e < a n (x) < M for all n <E N, x € R and 
sup ngN |z/„({a;})| < 1 — e for all Let us denote by X n the solution to ([3"T]) 

with a and v replaced by a n and v n , with respect to the same Brownian motion 
(and then on the same probability space) and assume that Xq converges in L 1 
to X Q . 

Assume moreover, that there exist two functions a and f in BV such that 
On > a and f Vn > / 
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and set v(dx) = —f'(dx)/(f(x) + f(x—)). Then X n converges uniformly in 
L X (R) to X on [0,T] for any T > 0, where X is the strong solution to X t = 
X + J* a{X s ) dB s + J R v{ dx)L% (X). 

Remark 15. This theorem endows the importance of the function / and allows us 
to recover the result of Proposition^ It has to be noted that one may construct 
a sequence (v n )n£N of measures in M satisfying the hypotheses of Theorem 1141 
converging weakly to a measure /x but for which the measure v given by this 
theorem differs from /i. Proposition [8] illustrates this points if we take for drift b 
a symmetric function with compact support on [—1, 1] with J_ b(x) — k/2, so 
that nb(x/n) dx converges to the Dirac measure (k/2)5q. Yet n~ 1 X tn 2 converges 
to the solution to Y t = B t + §^pj-L°(Y). 

J.-F. Le Gall also proved some extension of the Donsker theorem when a = 1 
in (l50|) . This Donsker theorem was used and improved by P. Etore in [2TJ, [28l [29] 
in order to approximate one-dimensional diffusion processes with discontinuous 
coefficients by random walks (See also Section 111.9.3)) . 

Recently R. Bass and Z. Q. Chen have shown in [5] strong existence and path- 
wise uniqueness of (|50| under the assumptions that i) v £ M, ii) a is bounded 
below and above by a positive constant iii) there exists a strictly increasing 
function / such that 

W(x) - a(y)\ 2 < \f(x) - f(y)\, Vx, y € M. (52) 

Let us note that if a G BV, then a satisfies (I5"2"j) with f(x) = 211(71100 \da\, 
so that the results of [9] are slightly more general than the one in [52] . 

The article [9] also provides us with a comparison principle and shows that 
there is no solution to (|50| when v({x}) > 1 for some x £ K (the non-existence 
of a solution was stated in |52j without proof). 

As the local time plays a central role in the theory of one-dimensional diffusion 
processes, SDE with local times are also considered in the study of diffusions 
with singular drift, in order to get some diffusions under weak hypotheses (see 
for example 0[lOl[3O], ...), Dirichlet processes [UES], ■•• In addition, they also 
appear as a tool to study the "critical cases" for results on strong/weak existence 
and uniqueness of solutions of SDE (see [23] [55] for example). 

Finally, let us note that this kind of process appears also naturally when one 
studies diffusion processes generated by differential operators with discontinuous 
coefficients: see Section 111 .91 

11.8.1. Spatial change of variable 

The symmetric Ito-Tanaka formula allows us to state a stability result of solu- 
tions of |50|) under a subclass of functions in BV. 

Proposition 17. Let F be a function whose derivative f belongs to BV and is 
bounded above and below by a positive constant constant. If X is the solution to 
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([50]) with a e BV and v G M, £/ien 

F(V t ) = F(V ) + / /(X>(X s )dB s 
Jo 



((/(*+) + /(*-)>( di) + /'( ds)) (53) 



// i* 1 is one-to-one, then Y — F(X) is also solution to some SDE of type ([50 
with a coefficient {fcr) o F^ 1 G BV and a measure fx G M wii/i 

/i( da;) = — — ^— g = ho F^ 1 where v{ Ax) - 



g{x+) + g(x— ) ' h(x) + h{x— ) 

for some function h. 

Proof. Formula (|53|) a direct consequence of the symmetric Ito-Tanaka formula 
(|32|) . After a one-to-one change of variable, if ip belongs to BV, then it is easily 
checked that J R ip'{ dx)k(x) = J R ip' ( dx)k(F~ 1 (x)) dx for any bounded, mea- 
surable function k. The measure \x is then computed with Lemma 2.1 in [52j 
which asserts the existence of h, the previous change of variable and the relation 
Lf 1 ^(X)^LnF(X)). □ 

As we saw it, the existence and uniqueness results on SDEs with local time 
relies on some application of this formula, to get rid of the local time. In [6] , 
M. Barlow used this kind of transform to construct counter-examples to the 
uniqueness of strong solutions to some SDEs. 

From a numerical point of view, this change of variable is important in order 
to simulate the process X, as shown in Section fl 1.9. 31 

In [67], Y. Ouknine made use of the Ito-Tanaka formula to study the distri- 
bution of solutions of SDEs of type (I50|) with a constant on R + and on R*_ and 
v = k5q. This class of SDEs remains stable when one uses for F in Proposition[T"7l 
functions of type F(x) = px + + qx~ . 



11.8.2. Random time change 

Random time change is also a tool to study SDE with local time, as shown first 
in [53] . In [ST] Chapter 5] , M. Martinez used a random time change based on the 
SBM instead of a Brownian motion to provide us with existence and uniqueness 
results of one-dimensional SDEs in the spirit of the work of H. Engelbert and 
J. Schmidt (see [21 or [46l Section 5.5.A]). 



11.9. Diffusions with discontinuous coefficients and diffusions on 
graphs 

We present now in this section some results about PDEs with discontinuous 
coefficients, and diffusion equation on graphs. This gives the general framework 
to consider both a large subclass of SDEs of type ([50]) and extensions of the 
Walsh's Brownian motion presented in Section 111.21 
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11.9.1. Infinitesimal generators and associated diffusion processes 

Let a and p be two measurable functions on R such that < A < a(x) < A and 
A < p{x) < A for all x £ R. Let also & be a measurable bounded 
Let (L,Dom(L)) be the differential operator 

L = Id! (°d!) + & cL? Dom W = {/eC(R;M)|L/eC(M;R)}. (54) 

Using scale functions and speed measures, or analytical results on the fundamen- 
tal solution of Jj — L, or Dirichlet forms, it can be proved that (L,Dom(L)) is 
the infinitesimal generator of a continuous, conservative, strong Markov process 
(X t ,t> 0;T t ,t> 0;P x ,xe R) (See [39l Ell [78] for example). 

The next proposition follows from an application of Corollary [2] and Theo- 
rem [TH (See SB ESI [54]). Note that a localization argument shall be used if b 
does not belong to L 1 (M;R) (see [751 Lemma II. 1.12] for example). 

Proposition 18. We assume in addition to the previous hypotheses that p and 
a belong to BV. Then for any starting point x, X is the unique strong solution 
to the SDE 

* 1 



X t =x+ I ^a{X s )p{X s )dB s + I -a'(X s )p(X s )d S 



o 



+ / b(X s )ds+ / V {dy)L\{X), (55) 



where B is a Brownian motion, L\{X) is the symmetric local time of X at the 
point y and 

V= Y1 ^7?ZV~7~3 ^" { Xn }«e Jcn = { * e R | a(x) ^ a(x-) } . 
a\x n -r i a\x n i 

The SDE (|55)) is a particular case of (|50)l encountered in Section fl 1.81 
The article [40j shows how to compute explicitly the density transition func- 
tion and the Green function of this process using spectral analysis. 



11.9.2. Transmission conditions and diffusions on a graph 

The solution of the parabolic PDE ^ = Lu with u(0,x) = f(x), where L is 
given by (|54p . has to be understood as a weak solution, that is as a function 
u{t,x) € C(0,T;L 2 (M)) n L 2 (0, T; H X (R)) such that for all ip € C°°([0, T];R), 

' d ^ t,X) u{t,x)p(x)- 1 dxdt- [ T [ b{x) dU ^ X) ^(t,x) dxdt 



dt Jq J b dx 



(i 



du(t,x)d1>(t,x) dxdt= _ f yM^o x ) p ( x )-i dx. 
ox OX /» 
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Now, if the coefficients a and p belong to BV and the points {xi}i e j at which 
a is discontinuous, then u is also solution to the transmission problem 

' d t u(t,x) = Lu(t,x) on R* X (R\{iEj } i6J ), 
^ u(0,a;) = f(x) on E, 
a{xi— )d x u(t, Xi— ) = a(a;i+)9 2; u(i, £j+) for i e J, 
,it(f, Xi— ) = u(t, for i € J, 



where tt G C^ 2 (R* + x (R \ { a* } ie 7 )) n C(R; x R) (See 09]). 
Remark 16. Note that we can get any arbitrary condition transmission at some 
point Xi without changing the diffusion coefficient on the intervals (xj_i, x{) and 
(xi,Xi+i). For this, it is sufficient to multiply the coefficients a by a pj and p 
by l/pj on (xj, Xj+i) for a well chosen sequence { ju,j }j eJ - 

Differential operators on graphs are natural expansions of one-dimensional 
differential operators with a transmission conditions at some given points. Let 
G be an oriented graph with a length t e > associated to each edge e. Thus 
each edge e may be seen as a segment of R and a differential operator L e — 
^a e (x)A + b e (x)V is associated to it. A differential operator (L,Dom(L)) may 
be constructed from the L e 's. Yet to properly define the operator (L,Dom(L)), 
one needs to specify the behavior of the functions in the domain Doin(L) at the 
vertices. The domain Dom(L) of L is composed of functions that are continuous 
on G, of class C 2 on each edge and such that at any vertex v, 

^ a v>e (v) = with a Vie € [0, 1], ^ a v,e = ^ ( 57 ) 

(J thy g 

where e ~ v means that v is an endpoint of the edge e and is the derivative 
of / in the direction of the edge e at the vertex v. 

The condition {57J prescribes the flux in each direction, and is clearly a 
transmission condition. The one-dimensional problem (|56[) can be seen as a 
diffusion equation on a graph whose vertices are the x^s and whose edges are 
the intervals (xi,Xi+i). 

Indeed this generalization goes farther, since one may derive an Ito formula 
for this diffusion process with a local time at each vertex [33]. Besides, ap- 
proximations by random walks is considered in [26j . The differential operator 
(L,Dom(L)) is the infinitesimal generator of a stochastic process X, which can 
be seen as the limit of a diffusion process moving in small tubes around the 
edges and reflected at the boundary: See [3S] . Diffusions on graphs are useful to 
describe the limiting behavior of perturbed Hamiltonian systems (see [33] and 
the subsequent works). Some of the Monte Carlo methods presented below in 
Section [T 1 . 9 .31 may be adapted to this case. On some particular graphs, it is also 
possible to compute explicitly the transition density function of X: see [BS] for 
example. 
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11.9.3. Numerical simulations 

The question of the simulation of the solution of some SDE with discontinuous 
diffusion coefficients is a problem which have been hardly treated (See however 
[45] [TBI [93]). Until recently, the case of processes generated by a divergence 
form operator of type hd Xi (aijd Xj ) have not been rigorously investigated, and 
the simulation of the process it generates in the multi-dimensional case remains 
an open problem (one knows from [50] that the diffusion may be approximated 
by a Markov chain, but the computation of the transition probabilities of this 
Markov chain is intractable in the general case). Due to many fields of applica- 
tions (geophysics, electro/magneto-encephalography, ecology, astrophysics, ...), 
it is however of practical importance, and the main problem is related to the 
understanding of the behavior of the diffusion process when it reaches a hyper- 
surface where a is discontinuous. 

In dimension one, Proposition [TH] together with Proposition [T7] allows us to 
describe exactly what happens to a diffusive particle that reaches a point where 
either p or a is discontinuous, if its trajectory is a realization of a trajectory of 
the process generated by the differential operator L given by (fM]) . Although one 
may think to use Corollary [2] to regularize the coefficients, this leads to unstable 
simulations (See [29] for numerical examples). Moreover, in practical situations, 
the jumps of the coefficients may be very high (for example, in a fissured porous 
media, the order of magnitude may be 1.000 or higher). In addition, it does not 
explain what happens to the diffusion. 

In dimension one, several schemes to simulate diffusion processes with discon- 
tinuous coefficients may be given, and we present here briefly several approaches. 
In many cases, we have assumed that the coefficients a and p belong to BV whose 
points of discontinuity have no cluster points, or that they can be approximated 
well by such coefficients. We also note that one may assume that b — by 
transforming the coefficients a and p into a(x) exp(— h(x)) and p(x) exp(h(x)) 
with h(x) = 2 Jq b(z) / 'a(z) p(z) dz. This is the Zvonkin transform, as we saw it 
in the proof of Proposition [8] 

We present several schemes. Once the behavior of the diffusion is understood, 
it becomes possible to provide other schemes or to mix them. 

A. If there is one discontinuity, one may easily adapt the method proposed by 
E. Hausenblas in [42] for the reflected diffusion processes. There, the simulation 
skips the small excursions when the process reaches the boundary, and the par- 
ticle jumps according to the entrance law of the diffusion process of excursions 
whose size is greater than a fixed parameter e. The adaptation can be done us- 
ing the decomposition of excursions measures given in [55] , and by choosing the 
sign of the excursion whose size is greater than e with an independent Bernoulli 
random variable. 

B. Using Proposition [171 one may find some one-to-one function F to get rid of 
the term with the local time in (|55|) . There, one obtains some SDE with a drift 
term and a diffusion term whose coefficients are discontinuous. It is then possible 
to apply results on the simulations of SDEs with discontinuous coefficients. This 
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is the strategy chosen by M. Martinez and D. Talay 61, 62J, where he uses the 
Euler scheme whose convergence is shown in [03] . In addition, the speed of 
convergence is computed. 

C. Still using Proposition [l7j after having approximated the coefficients by 
piecewise constant coefficients (after the drift have been removed), one may 
find a one-to-one function G that transform the solution X to ([55)) into the 
solution to the SDE 



Y t = G(X ) + B t + y: V^yPi^l- V^-)M*"-) L C(*") (n (58) 

where {x n } ne j is the set of points where a and p are discontinuous. Hence, 
around each point x n , Y behaves like a Skew Brownian motion. The idea 
is then the following: We construct a sequence of points {y k }keK such that 
{G(x n )} neJ C {y k }keK and t/ fe _i < y k < y k+1 . We then set 

e l = inf { t> e 1 - 1 1 Y t e {y k ~ 1 ,y k+1 } } 

if Ygi-i — y k . For this, we assume that 9° — and Ygo = Yq belong to {y k }keK- 
It is then possible to simulate exactly the Markov chain (8 l , Ygi ) . If Ygi = y k does 
not belong to {x n } ne j, then it corresponds to simulate the first exit time and 
position from [y k ~ 1 , y k ] for the Brownian motion, which can be done numerically 
with the help of analytical expressions of the density of the Brownian motion 
killed when it exit from some interval. If Ygi = y k = x n for n G J, then we 
have to assume that the {y k }keK are such that x n — y k ~ x = y k+1 — x n (this is 
why we use a finer grid {y k }ken)- Then, using the "symmetry" of the negative 
and positive excursions of the Skew Brownian motion, the first exit time and 
position for Y from [y ,y +1 ] can be deduced from the ones of the killed 
Brownian motion and the skewness parameter of Y around G{x n ). 

Let us note that it is also possible to simulate Yt for an arbitrary time 
T. For this, we decide at each step with a Bernoulli random variable whether 
{O 1 ^ 1 +8 l <T} or not. In the former case, it is also possible to simulate (O 1 , Ygi) 
given + l < T} when one knows Ygi-i). In the latter case, we are 

interested in simulating Yp given 8 l > T — which can also be done from 
the density of the Brownian motion killed when it exit from some interval. 

Once we have simulated the (ff 1 , ie*)i=o,i,..., it is easy to use G to get the 
distribution (r A T, X tA t), there r is the first exit time from some interval, and 
T > is an arbitrary time. 

This approach is presented in [53 |6T] . 

D. We have already noted in Section 111.9.31 that J.-F. Le Gall has proved a 
Donsker theorem of SDEs of type dZ t = dB t + J R v{ dx)Lf (Z). This means that 
he constructed a random walk (Sk)kefi such that n~ 1 S , ^ n 2 t j converges in distri- 
bution to Z t for any t > 0. The probability transitions of Sk can be explicitly 
deduced from the measure fi. 

In [27], P. Etore used the previous map G to reduce the simulation of X to 
this case, and then approximate the diffusion Y using the random walk. He also 
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computed the speed of convergence of this random walk. In some sense, this is 
also the simplification of the algorithm presented in C. where the {y k }keK are 
equally spaced (at size 1/n) and the exit times 8 1+1 — 8 % are replaced by their 
expectations, which are constant over all the position and equal to 1/n 2 . 

E. In [28 (see also [29), we constructed a bi-dimensional Markov chain (8 k , Z k )k£® 
from a realization of a trajectory of X generated by L by setting 

r k+1 =mf{t>T k \X t £G\ {X T u} } 
and 6 k+1 = 6 k +E[r fc+1 - r k | (X T h,X T h+i)], Z k+1 = X T k+i, 

where Q is an arbitrary grid whose distance between two points is not necessarily 
equal to a fixed parameter. 
Let us set 

z(t) = z k + o t ' 6 \ (Z k+1 -Z k ) 

for t € [8 k ,8 k+1 ] and K(t) = inf { k e N \ 8 k > t }. When the mesh of Q de- 
creases to 0, one may show that (8 Ki - t \ Z(t)) converges uniformly in t £ [0,T] 
in probability to (t,X t ) and the rate of convergence may be computed. This 
construction no longer uses the SDE (|55p . and the assumption that a and p 
belong to BV may be dropped. 

We then approximate X by simulating a Markov chain with the distribution 
of (Z k ,8 k )k=o,i.2,...- For this, the probability transition of (Z k )kefi is deduced 
from the scale function, since for 

gpr^i) - S{x) 

Vi{X) = — r — r for X € (Xi-i,X l+ i) 

S{x i+ i) - S(Xi-i) 

where Xi-\ < x.- L < Xi+i are three successive points on the grid G, then P[ Z k+1 — 
xi + i | Z k = %i] = Vi(Z k ). Let us note that Vi(x) is solution to Lvi(x) = on 
(xi-i, Xi+i) with «j(a;,_i) = 1 and Uj(a;,+i) = 1. The value of 8 k is then com- 
puted by solving Lui{x) = —Vi(x) on (a^-i, a^+i) with Ui(xi-i) = u^Xi+x) = 
andvi = Vi if Z k+1 = Xi+\ and Vi — 1 — Vi if Z k+1 = Xi-i. In facts, 8 k — Vi(Z k ). 

F. In [Hl[22], M. Decamps, A. De Schepper, M. Goovaerts and W. Shouten 
compute the densities of some process with discontinuous coefficients, in view 
of financial applications. For this, the Fokker-Planck equation is solved using 
perturbation formula of the density with respect to the density of the SBM: For 
the solution to the SDE 

dX t = a(X t ) AB t + b(X t ) dt + (2a - 1) dlf (X), 

where B is a Brownian motion, the write the density p(t, x, y) of X as 



P(t, x, y) = —I— exp [ f b(z) dz) pP(t, y ,y) 

\Jyo J 



xE yo 



exp(^-£ C (Xf)d S - 7 Lf (a; • 







(59) 
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where X 13 is the SBM(/3) with 

(g + l/2)tr(a;i) + (a - l/2)a(a:* ) 

; ~ (7(X* ) + (7«) 

p@(t,x,y) is its density, 

They then consider several "Skew models" — such as the Self Exciting Thresh- 
old (SET) Cox-Ingersoll-Ross, SET Vasicek, SET Libor market, ... — and use 
spectral decompositions to compute and/or approximate the density given in 

dSSD- 

They also consider the case of Skew Bessel processes. 



11.9.4- Parameter estimation 

In [TJini], O. Bardou and M. Martinez have constructed a scheme to estimate 
both the skewness parameter 7 and the position of the point x of the doubly 
reflected Skew Brownian motion dX t = dB t + jdLf(X) + dL^{X) - dL\(X). 
Their construction relies on the ergodicity of the underlying process, and its 
connections with PDEs with discontinuous coefficients. 



11.10. Multi- dimensional extensions 

The natural multi-dimensional extension of the Skew Brownian motion is that 
of a SDE involving the local time of the process in some hyper-surface. Yet it 
is difficult task to construct such a process. We then present briefly the method 
proposed by N. Portenko on generalized diffusion processes, where the semi- 
group of the process is constructed using a Volterra series. Afterwards, we give 
some references to other constructions, some of them relying on fixed point 
theorems on SDEs, and other on Dirichlet forms. 



11.10.1. Generalized diffusion processes 

In [69l [70l [71] , N. Portenko develops the concept of diffusions with singular drift, 
which we already used in Section [2l 

Let S be a surface in R d which separates R d into two regions: the interior 
region D and the exterior region R d \ D, N a vector field on S, q a continuous 
function from S to R. We assume that S satisfies both the interior and the 
exterior sphere property, which means that for any x S S, one can find some 
balls B and B' with B c D, B' c R d \ D, BnS = B'nS = {x}. 
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Let also a be a continuous function with values in the space of d x d-symmetric 
matrices, which is uniformly elliptic and bounded. A generalized diffusion is a 
diffusion process whose infinitesimal generator L is formally written 



d r^o d 



L = \H ^Wft^ + E^M^iC*)^. (60) 

ij — l 1 ^ i—1 1 

The idea is then to construct first the semi-group (Pt)t>o to L, it to show it 
generates a Feller process (X,¥ x ). The solution of the parabolic PDE 

^ = Lu(t, x) on M* + x R d with u(0, x) = <p(x) (61) 

is then u(t,x) = f Md p(t,x,y)ip(y) dy = E y [tp(X t )]. 

The problem is of course to interpret the meaning of the 5-function in the 
drift. One may see (|61| as a PDE with a transmission condition on S, as (|6ip 
may be written 

du{t,x) =Lu ^ x) onM rf\ 5j 

(1 + q{x))N+{x) ■ Vu(t, x) = {l- q(x))N-{x) ■ X7u{t, x) ( 62 ) 
u(t, •) is continuous on S, u(0, ■) = ip on K d , 

where N + and N_ are the inner and outer conormal to S, that is N + (x) = 
a(x)n(x) and N^(x) = —N + (x), where n(x) is the normal unit vector to S 
at x, which is directed inward. There, N±(x) ■ Vu(t,s) = lim e _,.o s~ (u(t, x ± 
eN±(x))-u(t,x)). 

As the 5-function imposes some flux condition on S, we can think to use a 
single layer potential, which is generated by a distribution of charges V(t, •) in 
S at time t which we choose to be 

V(t,x) = ^-(N + (x)-Vu(t,x+) + N-(x)-Vu(t,x-)), x e S. (63) 

Let p°(t,x,y) be the fundamental solution to L° = ^aijd x . Xj . The potential 
generated by V(t, x) on M. d is then equal to 

v(t,x)= / / p°{t - s,x,y)V(s,y)da y ds, 



o Js 



where a y is the Lebesgue measure on S, and according to classical results on 
potential theory, v is harmonic on K d \ S and continuous on R d (and then when 
crossing S). In addition, for x G S, 

N±(x)-Vv(t,x)= f f N+(x) .Vp°(t- s,x,y)V(s,y)da y dsTV(t,x). (64) 



As v(0, x) — 0, we then sought u to be the superposition of v and u (t, x) defined 

by 

u°(t,x) - / p a (t,x,y)<p(y)da y . (65) 
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This solution corresponds to the case q = and satisfies on S, N + (x)Vu°(t, x+) = 
N„(x)Vu°(t,x-). 

By combining the previous equation with Equations (|63[> . (|64|) . we obtain 
that in order that the flux condition in (f62|) is satisfied, the charge distribution 
V shall be solution to 



This equation may be solved using Volterra series, which means by constructing 
recursively 



and establishing the convergence of the series V(t, x) = J2n>o V (*) x ) ( tne 
advantage we get in Section [2] was that there was no need to consider such a 
series, since N+(x) ■ Vp°(t, 0,y) = when p°(t, x, y) is the one-dimensional heat 
kernel) . We have then constructed a charge distribution on S whose single layer 
potential u gives a solution to (|6"Tjl outside S. It is then easily obtained that u 
is solution to ([52]). 



We are now willing to construct the semi-group (Pt)t>o of L. For this, we set 
P t V(x) = J &d p°(t, x, y)ip(y) dy and 



where x € M. d , t > and V is the solution to (f66|) . This is a perturbation 
formula. It is then possible to prove that P s +ty>(x) = P s Pt(f{x) for any x £ K d , 
that Pt(p(x) > f>{ x ) f° r x S M d and that 



N+{x) ■ VP t f(x) = (1 -q(x))F(f,x) and A^_(x) • VP t ip{x) = (1 + q(x))V{t, x), 
so that (t,a;) i— > Ptf(x) is solution to ([62]) . 



As for the proof of Proposition [U if |<?(x)| < 1, then P t y; is non-negative 
when is non-negative. In addition P t l = 1, so that (Pt)t>o is a contraction 
semi-group for the uniform norm. 

Theorem 15 ([69, Theorem 1]). Under the previous hypotheses on S, a and q 
and if \q(x)\ < 1 for all x € S, The semi-group (Pt)t>o generates a conservative, 
continuous Markov process (X, ¥ x ,x £ W 1 ). 

Remark 17. If S is some hyper-plane in R" and a is the identity matrix, then 
similar computations also lead to the same result [69l Section 7]. 




(66) 
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One may wish to characterize the diffusion process constructed so by its 
diffusion and drift coefficients. It is standard that the diffusion coefficient is 
characterized by a(x)6 ■ 9 = lim^o E x [ ((X t - x) ■ 9) 2 ] for any 9 € M d . Using 
the semi- group and a continuous test function ip with compact support, one can 
show that 

/ tp(x)- [ {(y - x) ■ 9) 2 P{t 1 x, dy)dx >[ ip(x)a(x)6 ■ dx, 

JM d t JM d JM d 

where P(t, x, dy) is the probability transition of P. Still using the test function 
ip, one gets that is 

/ 1>(x)\ f {(y-x)-6)P(t,x,dy)dx—+ f iP(x)q(x)(N+(x) ■ 6) da x , 

which characterized a drift concentrated on S. 

In [70], N. Portenko decomposes the process X as a semi- martingale. 

Theorem 16 ( |701 Theorems 1 and 2]). There exists a continuous additive 
functional (Ct)t>o °f bounded variation such that M t = X t — X — Q is a con- 
tinuous square integrable martingale with {M) t — J^a(X s )ds and such that 
^ X [Q] = T t (q(x)NX(x)) and Q = f* l {XsSS} dC] for i = 1, . . . , d. 

Some of these results have been extended recently (see for example [72]), for 
example to take into account some elastic killing effect when the process passes 
through the surface S. 

11.10.2. SDE with local time in the multi- dimensional case 

There exists several works on SDEs with local time that may appears as gen- 
eralizations of the SBM. However, their constructions are rather technical in 
general, and the hypotheses tedious to write. This is why we give here just a 
few references. 

A. The article [51] provides a construction of some diffusion X associated to 
the Dirichlet form 

(u,v) £ Hj(K d ) ^ / a{x)Vu{x)Vv(x)(al G (x) + (1 - a)l M <,\ G (x))p(x) dx 

JM d 

where G is a domain of R d , a e (0, 1), a is uniformly elliptic, p > and of course, 
a and p have sufficient integrability condition. If a and p are "weakly differen- 
tiable" , then one can decompose X as some SDE with a local time associated 
to dG. 

B. Given a measure p which is singular with respect to the Lebesgue measure, Y. 
Oshima proved in [66] the weak existence and uniqueness of the d-dimensional 
process solution to 

d 

dX\ = J2 (<TiA X t) dB l + b i( X t) + nj{X t ) dM( + (3i(X t )ety ,i = l,...,d, 

3 = 1 
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where aa T = (a it j)f <j=1 , tt t = a/a 1A , a M > c > 0, bi = |Ejti^i a w' 
& = 2ai i Sj=2 9 Xj ctj,i, Pi = rij = Ti i = 0, -B is a ci-dimensional Brownian mo- 
tion, — J m £f /^(dxi), where (^ x )a;eR is a family of continuous non-decreasing 
processes such that If — f Q l{ X i=x} d^f an d 

/ £ x t f(x)d Xl = f f(Xl) ai>1 (X s )ds 
Jm Jo 

for all continuous function / with compact support, M is a family of continuous 
local martingales with {M\NP) t = 5 itj £% and (M\ B j ) = for i,j = 1, . . . , d. 
The coefficients a and r are assumed to be Lipschitz continuous. 

Y. Oshima also proves that this solution X is generated by the Dirichlet form 

(u,v) i— > — I a(x)Wu(x)Wv(x) dx + — j aij(x)Wu(x)Wv(x)i](dx), 

where ry is the measure r)( dx) = fi( dx\) dx2 ■ ■ ■ dx^, which is then the superpo- 
sition of two Dirichlet formio 

C. The previous construction was extended by S. Takanobu in [8TJ[82], where 
the assumptions that j3\ = and ai ( i > c > are replaced by the assumption 
that sup yeR d-i \/3d(y,x)fi({x})\ < 1 for all x £ K and Ylj=i a j,d > c > for 
some positive constant. 

D. In 77J, A.-S. Sznitman and S.R.S. Varadhan constructs the solution to the 
SDE 

N 

X t = x + B t + ^V k L° t {n k -X), 

k=0 

where the n^'s are unit vectors generating distinct hyperplanes H k , and V k are 
constant vectors with V k -n k = 0. This appears as a way to construct the solution 
of some N dimensional diffusion corresponding to N one-dimensional particles 
in interaction. 

E. In 95, 96, 97J, L. Zaitseva shows the strong existence and uniqueness of the 
SDE 

X t = x + Bt + / a{n s {X s )) dB Ls + f {(3v + a(n s (X s ))) dL s , (67) 
Jo Jo 

where S is a hyperplane of K d orthogonal to the unit vector v, its the projection 
from R d onto S, B a Brownian motion on M. , B a Brownian motion on S, 
G (— 1, 1), (L s ) s >o is the local time of v ■ X, a a measurable function from S 
to S and a a measurable function such that aa T takes its values in the space 
of symmetric matrices. The functions a and a are assumed to be bounded and 

2 The article 83 and subsequent works also contains some constructions of diffusion pro- 
cesses associated to the superposition of Dirichlet forms defined on hypersurfaces, but without 
giving some SDE solved by them. 
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Lipschitz continuous. This diffusion is shown in 97J to be a generalized diffusion 
in the sense of N. Portenko [7T] with drift (fiv + ao tts)Ss and a diffusion term 
Id + a o ttsSs- 

11.11. Domains of applications 

Being related to diffusion processes with discontinuous coefficients and to diffu- 
sions on graphs, the Skew Brownian motion has potential applications in many 
fields of physics: in astrophysics [98], in geophysics and study of heterogeneous 
media [35] [SO] [56] Ell EH] [92] , in perturbed Hamiltonian systems [37] , study of 
hysteresis phenomena [32j , in biology (some modelling of nerve impulses involves 
parabolic PDEs as in Section fl 1.9. 21 [51] ). ecology (the dynamic of populations 
moving between different refuges in studied in |15|). in finance and actuarial 
sciences [HI [20j [21] , in optimization (see below)... 

In all these applications, the properties of the Skew Brownian motion may 
be used in order to solves the related problems by Monte Carlo methods, or to 
use its properties to solve such a problem analytically. 

Applications to optimal problems 

The article [53] looks for a a drift b > with L b(x) dx < +oo that minimizes 
Ex t, where r = inf { t > | X T = 1 } and X is the solution to dX t — dB t + 
b(X t ) dt with a reflecting boundary at 0. Indeed, the optimal drift is a measure 
that can be computed explicitly, and that is of type b(dx) = b(x)dx + kS Xo (dx). 
Hence it involves a process of Skew Brownian motion type. 

The constructions in [5] and in [T7] are also related to methods for solving 
resources' allocations problems (see [5U] for example). Recently, the variably 
skew Brownian motion presented in Section [11. 61 has also been used for solving 
the Skorokhod embedding problem [18] . 

Acknowledgment. The author is grateful to Pierre Etore for interesting dis- 
cussions about the Skew Brownian motion, especially about the approach by 
N. Portenko. The author also wishes to thank the referees who helped him to 
improve the quality of this article. 

References 

[1] Bardou, O. (2005). Controle dynamique des erreurs de simulation et 
d'estimation de processus de diffusion. Ph.D. thesis, Universite de Nice 
and INRIA Sophia- Antipolis. 

[2] Barlow, M., Burdzy, K., Kaspi, H. and Mandelbaum, A. (2000). Variably 
skewed Brownian motion. Electron. Comm. Probab. 5, 57-66. MR1752008 

[3] Barlow, M., Burdzy, K., Kaspi, H. and Mandelbaum, A. (2001). Coales- 
cence of skew Brownian motions. In Seminaire de Probabilites, XXXV. 
Lecture Notes in Math., Vol. 1755. Springer, 202-205. [MR1837288I 



A. Lejay/The skew Brownian motion 



461 



[4] Barlow, M., Pitman, J. and Yor, M. (1989a). On Walsh's Brownian mo- 
tions. In Seminaire de Probabilites, XXIII. Lecture Notes in Math., Vol. 
1372. Springer, 275-293. IMR.1022917I 

[5] Barlow, M., Pitman, J. and Yor, M. (1989b). Une extension multidimen- 
sionnelle de la loi de l'arc sinus. In Seminaire de Probabilites, XXIII. 
Lecture Notes in Math., Vol. 1372. Springer, 294-314. [MR1022918I 

[6] Barlow, M. T. (1988). Skew Brownian motion and a one-dimensional 
stochastic differential equation. Stochastics 25, 1, 1-2. MR1008231 

[7] Barlow, M. T., Emery, M., Knight, F. B., Song, S. and Yor, M. (1998). 
Autour d'un theoreme de Tsirelson sur des nitrations browniennes et 
non browniennes. In Seminaire de Probabilites, XXXII. Lecture Notes 
in Math., Vol. 1686. Springer, 264-305. 

[8] Bass, R. and Chen, Z.-Q. (2002). Stochastic differential equations for 
dirichlet processes. Probab. Theory Related Fields 121, 3, 422-446. 
<doi:10 . 1007/s004400100151>. MR1867429I 

[9] Bass, R. and Chen, Z.-Q. (2005). One-dimensional stochastic differential 
equations with singular and degenerate coefficients. Sankhyd 67, 19-45. 
IMR2203887I 

[10] Bass, R. F. and Chen, Z.-Q. (2003). Brownian motion with singular drift. 

Ann. Probab. 31, 2, 791-817. IMR1964949I 
[11] Blumenthal, R. (1992). Excursions of Markov Processes. Probability and 

Its Applications. Birkhauser. MR113846T1 
[12] Breiman, L. (1981). Probability. Addison- Wesley. 

[13] Burdzy, K. and Chen, Z.-Q. (2001). Local time flow related to skew Brow- 
nian motion. Ann. Probab. 29, 4, 1693-1715. MR1880238 

[14] Burdzy, K. and Kaspi, H. (2004). Lenses in skew Brownian flow. Ann. 
Probab. 32, 4, 3085-3115. MR2094439] 

[15] Cantrell, R. and Cosner, C. (1999). Diffusion models for population dy- 
namics incorporating individual behavior at boundaries: Applications to 
refuge design. Theoritical Population Biology 55, 2, 189-207. 

[16] Chan, K. S. and Stramer, O. (1998). Weak consistency of the Euler method 
for numerically solving stochastic differential equations with discontinuous 
coefficients. Stochastic Process. Appl. 76, 1, 33-44. MR16380H] 

[17] Cherny, A., Shiryaev, A. and Yor, M. (2004). Limit behavior of the 
"horizontal-vertial" random walk and some extension of the Donsker- 
Prokhorov invariance principle. Theory Probab. Appl. 47, 3, 377-394. 

[18] Cox, A.M.G. and Hobson, D.G. (2005). A unifying class of Skorokhod 
embeddings: connecting the Azema-Yor and Vallois embeddings. Preprint, 
< |arxiv:math.PR/0506040^ >.[ MR2213778l 

[19] Decamps, M., de Schepper, A. and Goovaerts, M. (2004). Applications of 
J-perturbation to the pricing of derivative securities. Phys. A 342, 3-4, 
677-692. IMR2105796I 

[20] Decamps, M., de Schepper, A., Goovaerts, M. and Schoutens, W. (2005). 
A note on some new perpetuities. Scand. Actuar. J. 2005, 4, 261-270. 
<DOl:10 . 1080/03461230510009772>. 

[21] Decamps, M., Goovaerts, M. and Schoutens, W. (2006). Asymetric 



A. Lejay/The skew Brownian motion 



462 



skew Bessel processes and related processes with applications in finance. 
J. Comput. Appl. Math. 186, 1, 130-147. 1MR2190302I 
Decamps, M., Goovaerts, M. and Schoutens, W. (2006). Self Exciting 
Threshold Interest Rates Model. Int. J. Theor. Appl. Finance 9, 7, 1093- 
1122, <DOl:10 . 1142/S0219024906003937>. 

Desbois, J. (2002). Occupation times distribution for Brownian motion 
on graphs. J. Phys. A: Math. Gen. 35, L673-L678. [MR1947234I 
Engelbert, H. J. and Schmidt, W. (1985). On one-dimensional stochas- 
tic differential equations with generalized drift. In Stochastic differential 
systems (Marseille-Luminy, 1984)- Lecture Notes in Control and Inform. 
Sci., Vol. 69. Springer, 143-155. [MR0798317I 

Engelbert, H.-J. and Wolf, J. (1999). Strong Markov local Dirichlet pro- 
cesses and stochastic differential equations. Theory Probab. Appl. 43, 2, 
189-202. MR 1679006I 

Enriquez, N. and Kifer, Y. (2001). Markov chains on graphs and Brownian 
motion. J. Theoret. Probab. 14, 2, 495-510. IMR1838739I 
Etore, P. (2006). On random walk simulation of one-dimensional diffusion 
processes with discontinuous coefficients. Electron. J. Probab. 11, 9, 249- 
275. MR 2217816I 

Etore, P. and Lejay, A. (2006). A Donsker theorem to simulate one- 
dimensional diffusion processes with measurable coefficients. Preprint, 
Institut Elie Cartan, Nancy, France. 

Etore, P. (2006). Approximation de processus de diffusion a coefficients 
discontinus en dimension un et applications a la simulation. Ph.D. thesis, 
Universite Nancy 1. 

Flandoli, F., Russo, F. and Wolf, J. (2003). Some SDEs wit h distribution al 

drift. I. General calculus. Osaka J. Math. 40, 2, 493-542. [MR1988703I 

Folland, G. B. (1995). Introduction to partial differential equations, 2 ed. 

Princeton University Press, Princeton, NJ. MR1357411 

Freidlin, M., Mayergoyz, I.D. and Pfeifer, R. (2000). Noise in hysteretic 

systems and stochastic processes on graphs. Phys. Rev. E 61, 2, 1850- 

1855. 

Freidlin, M. and Sheu, S.-J. (2000). Diffusion processes on graphs: stochas- 
tic differential equations, large deviation principle. Probab. Theory Related 
Fields 116, 2, 181-220. MR1 7437691 

Freidlin, M. and Weber, M. (1998). Random perturbations of nonlinear 
oscillators. Ann. Probab. 26, 3, 925-967. [MR1634409I 
Freidlin, M. and Wentzell, A. (1993). Diffusion processes on graphs and 
the averaging principle. Ann. Probab. 21, 4, 2215-2245. IMR12453081 
Freidlin, M. and Wentzell, A. (1994a). Necessary and sufficient conditions 
for weak convergence of one-dimensional Markov processes. In Festschrift 
dedicated to 70 th Birthday of Professor E.B. Dynkin, M. Freidlin, Ed. 
Birkhauser, 95-109. [MRi3I1713| 

Freidlin, M. I. and Wentzell, A. D. (1994b). Random perturbations of 
Hamiltonian systems. Mem. Amer. Math. Soc. 109, 523. MR 12012691 
Friedman, A. (1983). Partial Differential Equations of Parabolic Type. 



A. Lejay/The skew Brownian motion 



463 



Robert E. Krieger Publishing Company. MR0454266 
[39] Fukushima, M., Oshima, Y. and Takeda, M. (1994). Dirichlet Forms and 

Symmetric Markov Process. De Gruyter. MR1303354 
[40] Gavcau, B., Okada, M. and Okada, T. (1987). Second order differential 

operators and Dirichlet integrals with singular coefficients. I: Functional 

calculus of one-dimensional operators. Tdhoku Math. J. (2) 39, 4, 465-504. 

IMR09174631 

[41] Harrison, J. and Shepp, L. (1981). On skew Brownian motion. Ann. 
Probab. 9, 2, 309-313. IMR0606993I 

[42] Hausenblas, E. (1999). A numerical scheme using ltd excursions for sim- 
ulating local time resp. stochastic differential equations with reflection. 
Osaka J. Math. 36, 1, 105-137. MR99m:60087. IMR16707541 

[43] Ikeda, N. and Watanabe, S. (1989). Stochastic differential equations 
and diffusion processes 2 ed. North-Holland Mathematical Library 24. 
IMR10112521 

[44] Ito, K. and McKean, H. (1974). Diffusion and Their Sample Paths, 2 ed. 
Springer- Verlag. 

[45] Janssen, R. (1984). Difference methods for stochastic differential equations 
with discontinuous coefficients. Stochastics 13, 3, 199-212. MR0762816 

[46] Karatzas, I. and Shreve, S. (1991). Brownian Motion and Stochastic Cal- 
culus, 2 ed. Springer- Verlag. IMR1121940I 

[47] Kurtz, T. and Protter, P. (1996). Weak convergence of stochastic in- 
tegrals and differential equations. In Probabilistic Models for Nonlinear 
Partial Differential Equations, Montecatini Terme, 1995, D. Talay and 
L. Tubaro, Eds. Lecture Notes in Math., Vol. 1627. Springer- Verlag, 1- 
41. MR1431298I 

[48] Ladyzenskaja, O., Solonnikov, V. and Ural'ceva, N. (1968). Linear and 
Quasilinear Equations of Parabolic Type. American Mathematical Society. 

[49] Ladyzenskaja, O. A., Rivkind, V. J. and Ural'ceva, N. N. (1966). Solv- 
ability of diffraction problems in the classical sense. Trudy Mat. Inst. 
Steklov. 92, 116-146. MR0211050I 

[50] Lang, R. (1995). Effective conductivity and Skew Brownian motion. 
J. Statist. Phys. 80, 125-145. [MR1340556I 

[51] Lc Gall, J.-F. (1982). Applications du temps local aux equations 
diffcrcnticllcs stochastiques unidimensionnellcs. In Seminaire de Proba- 
bilites XVII. Lecture Notes in Math., Vol. 986. Springer- Verlag, 15-31. 

[52] Le Gall, J.-F. (1985). One-dimensional stochastic differential equations 
involving the local times of the unknown process. In Stochastic Analysis 
and Applications. Lecture Notes in Math., Vol. 1095. Springer- Verlag, 
51-82. MR0 777514I 

[53] Lee, S. (1994). Optimal drift on [0,1]. Trans. Amer. Math. Soc. 346, 1, 

159-175. 1MR1254190I 
[54] Lejay, A. (2000). Methodes probabilistes pour Phomogeneisation des 

operateurs sous forme-divergence : cas lineaires et semi-lineaires. Ph.D. 

thesis, Universite de Provence, Marseille, France. 
[55] Lejay, A. (2001). On the decomposition of excursions measures of processes 



A. Lejay/The skew Brownian motion 



464 



whose generators have diffusion coefficients discontinuous at one point. 
Markov Process. Related Fields 8, 1, 117-126. 1MR1897608I 
Lejay, A. (2003). Simulating a diffusion on a graph, application to reservoir 
engineering. Monte Carlo Methods Appl. 9, 3, 241-256. [MR2009371I 
Lejay, A. (2004). Monte Carlo methods for fissured porous media: a grid- 
less approach. Monte Carlo Methods Appl. 10, 3-4, 385-392. Conference 
proceeding of IV IMACS Seminar on Monte Carlo Methods, . IMR21 05066 
Lejay, A. and Martinez, M. (2006). A scheme for simulating one- 
dimensional diffusion processes with discontinuous coefficients. Ann. 
Appl. Probab. 16, 1, 107-139. <DOi:10 . 1214/105051605000000656>. 
MR2209338J 

Ma, Z. and Rockner, M. (1991). Introduction to the Theory of (Non- 
Symmetric) Dirichlet Forms. Universitext. Springer- Verlag. MRU 06853 
Mandelbaum, A., Shepp, L. A. and Vanderbei, R. J. (1990). Optimal 
switching between a pair of Brownian motions. Ann. Probab. 18, 3, 1010- 
1033. IMRT062057I 

Martinez, M. (2004). Interpretations probabilistes d'operateurs sous forme 
divergence et analyse de methodes numeriques associees. Ph.D. thesis, 
Universite de Provence / INRIA Sophia- Antipolis. 

Martinez, M. and Talay, D. (2006). Discretisation d'equations 
differentielles stochastiques unidimensionnelles a generateur sous forme di- 
vergence avec coefficient discontinu. C. R. Math. Acad. Sci. Paris 342, 1, 
51-56. [MR2193396I 

Nakao, S. (1972). On the pathwise uniqueness of solutions of one- 
dimensional stochastic differential equations. Osaka J. Math. 9, 513-518. 
IMR03268401 

Nicaise, S. (1985). Some results on spectral theory over networks, applied 
to nerve impulse transmission. In Orthogonal polynomials and applications 
(Bar-le-Duc, 1984). Lecture Notes in Math., Vol. 1171. Springer, 532-541. 
MR0 8390241 

Okada, T. (1993). Asymptotic behavior of skew conditional heat kernels 
on graph networks. Canad. J. Math. 45, 4, 863-878. MR 12276641 
Oshima, Y. (1982). Some singular diffusion processes and their associated 
stochastic differential equations. Z. Wahrsch. Verw. Gebiete 59, 2, 249- 
276. IMR065"06"l6l 

Ouknine, Y. (1990). "Skew-Brownian motion" and derived processes. The- 
ory Probab. Appl. 35, 1, 163-169. IMR1050069I 

Peskir, G. (2005). A change-of-variable formula with local time on curves. 
J. Theoret. Probab. 18, 3, 499-535. <DOi:10 . 1007/sl0959-005-3517-6>. 
IMR21676401 

Portenko, N. (1979a). Diffusion processes with a generalized drift coeffi- 
cient. Theory Probab. Appl. 24, 1, 62-78. [MR0522237I 
Portenko, N. (1979b). Stochastic differential equations with generalized 
drift vector. Theory Probab. Appl. 24, 2, 338-353. [MR0532446I 
Portenko, N. I. (1990). Generalized diffusion processes. Translations of 
Mathematical Monographs, Vol. 83. American Mathematical Society, 



A. Lejay/The skew Brownian motion 



465 



Providence, RI. MR1 1046601 

[72] Portenko, N. I. (2001). A probabilistic representation for the solution to 
one problem of mathematical physics. Ukrainian Math. J. 52 (2000), 9, 
1457-1469. IMR18169401 

[73] Ramirez, J.M., Thomann, E.A., Waymire, E.C., Haggerty R. and Wood B. 
(2005). A generalization Taylor- Aris formula and Skew Diffusion Multi- 
scale Model. Simul. 5, 3, 786-801, <DOI:10 . 1137/050642770>. 

[74] Rogers, L. and Williams, D. (2000). ltd Calculus, 2 ed. Diffusions, Markov 
Processes, and Martingales, Vol. 2. Cambridge University Press. 

[75] Rosenkrantz, W. A. (1974/75). Limit theorems for solutions to a class 
of stochastic differential equations. Indiana Univ. Math. J. 24, 613-625. 
IMR03681431 

[76] Rozkosz, A. (1996). Weak convergence of diffusions corresponding to 
divergence form operator. Stochastics Stochastics Rep. 57, 129-157. 
IMR140795T1 

[77] Sznitman, A.-S. and Varadhan, S. R. S. (1986). A multidimensional pro- 
cess involving local time. Probab. Theory Relat. Fields 71, 4, 553-579. 

IMR08332691 

[78] Stroock, D. (1988). Diffusion semigroups corresponding to uniformly ellip- 
tic divergence form operator. In Seminaire de Probabilites XXII. Lecture 
Notes in Math., Vol. 1321. Springer- Verlag, 316-347. IMR0960535I 

[79] Stroock, D. and Varadhan, S. (1979). Multidimensional diffusion pro- 
cesses. Grundlehrcn der Mathematischcn Wissenschaften, Vol. 233. 
Springer- Verlag, Berlin. MR0 532498I 

[80] Stroock, D. and Zheng, W. (1997). Markov chain approximations to sym- 
metric diffusions. Ann. Inst. H. Poincare Probab. Statist. 33, 5, 619-649. 
IMR14735681 

[81] Takanobu, S. (1986). On the uniqueness of solutions of stochastic dif- 
ferential equations with singular drifts. Publ. Res. Inst. Math. Sci.11, 5, 
813-848. MR0866659 

[82] Takanobu, S. (1987). On the existence of solutions of stochastic differential 
equations with singular drifts. Probab. Theory Related Fields 74, 2, 295- 
315. MR0871257I 

[83] Tomisaki, M. (1980). Superposition of diffusion processes. J. Math. Soc. 

Japan 32, 4, 671-696. IMR0589106I 
[84] Trutnau, G. (2005). Multidimensional skew reflected diffusions. Stochastic 

analysis: classical and quantum, 228-244, World Sci. Publ., Hackensack, 

NJ. MR2233163I 

[85] Tsirelson, B. (1997). Triple points: from non-Brownian nitrations to har- 
monic measures. Geom. Fund. Anal. 7, 6, 1096-1142. MR1 4877551 

[86] Vuolle-Apiala, J. (1996). Skew Brownian motion-type of extensions. 
J. Theor. Probab. 9, 4, 853-861. [MR1419866I 

[87] Walsh, J. (1978). A diffusion with a discontinuous local time. In Temps 
locaux. Asterisque. Societe Mathematique de France, 37-45. MR0509476 

[88] Walsh, J. (1981). Optional increasing paths. In Two-index random pro- 
cesses (Pans, 1980). Lecture Notes in Math., Vol. 863. Springer, 172-201. 



A. Lejay/The skew Brownian motion 



466 



MR0630313 

[89] Watanabe, S. (1984). Excursion point processes and diffusions. In Proceed- 
ings of the International Congress of Mathematicians (Warsaw, 1983), 
Vol. 1, 2, 1117-1124, PWN, Warsaw. [MR0804763I 

[90] Watanabe, S. (1995). Generalized arc-sine laws for one-dimensional dif- 
fusion processes and random walks. In Stochastic analysis (Ithaca, NY, 
1993). Proc. Sympos. Pure Math., Vol. 57. Amer. Math. Soc, 157-172. 
IMR13354701 

[91] Weinryb, S. (1982). Etude d'une equation differentielle stochastique avec 

temps local. In Seminaire de Probabilites XVII. Lecture Notes in Math., 

Vol. 986. Springer- Verlag, 72-77. IMR0770397I 
[92] Weinryb, S. (1984). Homogeneisation pour des processus associes a des 

frontieres permeables. Ann. Inst. H. Poincare Probab. Statist. 20, 4, 373- 

407. IMR0771896I 

[93] Yan, L. (2002). The Euler scheme with irregular coefficients. Ann. 

Probab. 30, 3, 1172-1194. MR1920104I 
[94] Yor, M. (1997). Some aspects of Brownian motion. Part II: Some recent 

martingale problems. Lectures in Mathematics ETH Zurich. Birkhauser 

Verlag. MR1442263J 
[95] Zaitseva, L. (2000). On a probabilistic approach to the construction of the 

generalized diffusion processes. Theory Stoch. Process. 6, 1-2, 141-146. 

IMR20077321 

[96] Zaitseva, L. (2005). On stochastic continuity of generalized diffusion pro- 
cesses constructed as the strong solution to an SDE. Theory Stoch. Pro- 
cess. 11(27), 1-2, 125-135. < |arxiv:math.PR/0609305| >. |MR2QQ773l 

[97] Zaitseva, L. (2005). On the Markov property of strong solutions to SDE 
with generalized coefficients. Theory Stoch. Process. 11(27), 3-4, 140-146. 
< |arxiv:math.PR/0609307f >. 

[98] Zhang, M. (2000). Calculation of diffusive shock acceleration of charged 
particles by skew Brownian motion. Astrophys. J. 541, 428-435. 

[99] Zhikov, V., Kozlov, S. and Oleinik, O. (1981). G-convergence of parabolic 
operators. Russian Math. Survey 36, 1, 9-60. IMR0608940I 
[100] Zhikov, V., Kozlov, S., Oleinik, O. and T'en Ngoan, K. (1979). Averaging 
and G-convergence of differential operators. Russian Math. Survey 34, 5, 
69-147. MR0562800 



